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1 .   INTRODUCTION 

In  a  recent  article  Foutz  (1980)  introduced  a  new  test  for  goodness-of- 

fit,  to  be  called  the  F   test  in  the  sequel.   Although  the  test  was  proposed 

for  fitting  a  continuous  p-variate  distribution,  it  applies  equally  well  to 

univariate  problems.   The  null  distribution  of  the  test  statistic  was  shown 

to  be  distribution  free  as  well  as  being  independent  of  the  number  of  vari- 

ates  p.   Foutz  obtained  an  integral  representation  for  the  null  CDF  of  F  : 

n 

explicit  expressions  for  this  CDF  were  given  for  sample  sizes  2  and  3. 

Closed  form  solutions  for  the  CDF  for  larger  sample  sizes  are  quite  hard 

to  derive  and  Foutz  has  provided  a  large  sample  normal  approximation  to  the 

null  distribution  of  F  . 

n 

In  a  preliminary  comparison  with  ten  replications  of  50  simulated  samples 
from  (I)  a  mixture  of  uniform  distributions  and  @.i)  a  standard  normal  distri- 
bution, Foutz  found  that  the  F   test  outperformed  both  the  Chi  squared  test 

n 

and  the  Kolmogorov-Smirnov  (K-S)  test. 

In  this  paper  we  present  the  results  of  an  extensive  investigation  to 

compare  the  three  goodness-of-f it  tests,  the  Chi  squared  test,  the  K-S  test 

and  the  F  test.   Members  from  three  families  of  distributions,  viz.,  the 
n 

family  of  asymmetric  stable  distributions,  mixtures  of  normals,  and  the 
Pearson  family  have  been  selected  to  represent  the  true  underlying  distri- 
bution of  the  samples.  The  goodness-of-f it  tests  are  applied  to  test  the 
hypothesis  that  the  samples  are  from  a  standard  normal  distribution.   The 
measure  of  comparison  used  in  the  study  is  the  empirical  power,  based  on 
5000  replications  of  each  of  the  tests. 

The  Chi  squared  and  Kolmogorov-Smirnov  statistics  were  computed  using 
the  IMSL  routines  GFIT  and  NKS1.   The  number  of  cells  used  in  the  Chi 


International  Mathematical  and  Statistical  Libraries,  7500  Bellaire 
Boulevard,  Houston,  TX  77036 


squared  computation  was  6,  8,  and  8  for  sample  sizes  of  20,  30,  and  50, 
respectively. 

A  brief  discussion  of  the  F  test  is  given  in  Section  2  and  a  descrip- 
tion of  the  simulation  is  in  Section  3.   The  results  of  the  simulation  are 
presented  in  Section  4.   FORTRAN  codes  used  for  the  simulation  and  detailed 
tables  of  simulation  results  are  in  Appendices  I  and  II. 


2.   F  TEST 
n 

The  F   test  is  based  on  a  comparison  of  a  continuous  empirical  distri- 
n 

bution  function  (CEDF)  with  the  hypothesized  CDF.   The  CEDF  is  obtained  by 
"spreading"  the  total  mass  over  "statistically  equivalent  blocks"  generated 
by  the  sample.   As  shown  in  Anderson  (1966)  and  Foutz  (1980),  given  the 
order  statistics  of  a  random  sample  of  size  (n-1) ,  n  statistically  equiva- 
lent blocks  that  partition  the  sample  space  can  be  constructed  in  many 
different  ways  by  choosing  what  are  called  cutting  functions.   An  intuitively 
appealing  set  of  blocks,  which  is  the  one  used  in  this  study  is  obtained  by 
choosing  the  identity  functions  for  the  cutting  functions.   In  this  case, 
the  n  statistically  equivalent  blocks  are  B  ■  (-<»,  x,^],   B  =  (x.  .,x.  .],, 

B  =  (x,   1N)»)  where  x/1N,  x/ox,...,  x.   ,N  are  the  order  statistics  of 
n     (n-1)  (.1)    U)       (n-1) 

a  sample  of  size  n-1.   The  CEDF  is  constructed  by  spreading  a  mass  —  cont- 
inuously and  in  the  same  proportion  as  the  hypothesized  CDF  over  each 

block.   If  H  is  the  hypothesized  CDF  and  H  the  CEDF,  the  test  statistic 
o  n 

F  is  defined  as 
n 

F  =  SuplH  (x)  -  H  (x) I  (1) 

n    xr  '  n       o    ' 

Let  D.,  i  =  l,2,...,n  be  the  probability  contents  of  the  blocks  B. 

under  the  null  CDF  H  ,  i.e.,  D.  =  P[x  e  B.lH  ].   A  computationally  convenient 

o         1  l   o 

form  for  F  can  be  shown  to  be 
n 

n         1 
F  =  I     max(0,  -  -  D.)  .  (2) 

n   i-1        n    X 

Foutz  has  provided  the  null  distribution  of  F  in  integral  form  and  derived 

n 

closed  form  solutions  for  n  =  3,4.   Simplifying  the  expressions  given  by 
Foutz  for  n  =  3,4  yields 


P(F3   <   x)    = 


6x 


1  -   3(| 


-   x)' 


0    <   x  <  - 
|  <   x   <   2/3 
x  >   2/3 


(3) 


and 


/    20x3 


P(F     <  x)   = 

4    - 


3  2        9  1 

-20x     +   18x     -  fx    +  tt 
4  16 

3  3 

1   -   4(f  -  x)J 


I    * 


0  <  x  <  -r 

-        -  4 

1  1 
4±X<  2 


We  have  also  obtained  the  expression  for  the  null  CDF  for  n  =  5, 

70x4 

-105x4  +  80x3  -  12x2  +  |yx  -  yIj 

P(F,  <  x)  =(  .-  4   Qn  3   228  2   176    31 
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(5) 


As  is  evident  the  exact  distribution  is  quite  difficult  to  obtain  for 
higher  sample  sizes.   A  large  sample  normal  approximation,  due  to  Foutz  is 
given  by 


lim  P 
n-*» 


F  <  x 
n  - 


=  $ 


n(x  -  e  ) 

-1     -2  h 
(2e  X  -  5e  ) 2 


(6) 


where  $  is  the  standard  normal  CDF.   For  n-1  ■  20,30,50  we  used  this  approxi- 
mation to  test  the  hypothesis  that  a  simulated  sample  from  U[0,1],  a  uniform 
distribution  on  [0,1],  is  in  fact  from  that  distribution.   In  80,000  repli- 
cations, the  observed  significance  level  (number  of  hypothesis  rejections/ 
80,000)  was  consistently  smaller  than  the  nominal  value  as  can  be  seen  from 


TABLE  1 

EMPIRICAL  SIGNIFICANCE  LEVEL  OF 

IT  USING  ASYMPTOTIC  AP] 

(80,000  REPLICATIONS) 


FOUTZ  Fn  TEST  USING  ASYMPTOTIC  APPROXIMATION 


Sample  Size  20  30  50 

Significance  Level 
.10 
.05 
.01 


0757 

.0800 

.0859 

0372 

.0399 

.0428 

0082 

.0083 

.0093 

Table  1.   We  therefore  constructed  a  Monte  Carlo  CDF  of  F  for  n-1  =  2,3,4, 

n 

20,30,50  based  on  25,000  computer  generated  F  values;  these  represent  values 

n 

of  the  F  statistic  for  testing  the  hypothesis  that  a  set  of  samples  from 

U[0,1]  is  in  fact  from  U[0,1].   A  comparison  of  the  Monte  Carlo  CDF  with 

the  exact  CDF  for  n-1  =  2,3,4  is  provided  in  Table  2.   It  can  be  seen  that 

the  Monte  Carlo  CDF  provides  a  reasonable  approximation  even  for  small  n. 

The  power  properties  of  the  F  test  detailed  in  this  paper  are  based  on 

the  Monte  Carlo  CDF  of  F  .   Critical  values  obtained  from  the  Monte  Carlo 

n 

simulation  for  significance  levels  .01,  .05,  .1  and  n-1  =  20,30,50  are  in 
Table  3.  In  Table  4  we  present  the  observed  significance  level  in  225,000 
replications  when  testing  if  a  set  of  samples  from  U[0,1]  is  in  fact  from 
U[0,1]. 
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TABLE  2 
MONTE  CARLO  SIMULATION  VS   EXACT  VALUES  OF  CDF 


.40 
.45 
.50 
.55 
.60 
.65 


n=3 
MC         EXACT 


,7883 
.8604 
,9182 
,9598 
,9864 
,9989 


.7867 
.8592 
.9167 
.9592 
.9867 
.9991 


n=4 
MC         EXACT 


7624 
8725 
9396 
9694 
9860 
9952 


.7625 
.8725 
.9375 
.9680 
.9865 
.9960 


n=5 
MC         EXACT 


7607 
,8671 
,9405 
,9768 
,9915 

9975 


.7600 
.8693 
.9405 
.9778 
.9920 
.9975 
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TABLE  3 


CRITICAL  VALUES  FOR  Fn  TEST 
OBTAINED  BY  MONTE  CARLO  SIMULATION 


Sample  Size  20  30  50 

Significance  Level 
.10 
.05 
.01 


42714 

.41903 

.40816 

44865 

.43553 

.42116 

48659 

.46579 

.44487 
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TABLE  4 

EMPIRICAL  SIGNIFICANCE  LEVEL  OF 

1ST  USING  MONTE  CARLO  Al 

(225,000  REPLICATIONS) 


FOUTZ  Fn  TEST  USING  MONTE  CARLO  APPROXIMATION 


Sample  Size  20  30  50 

Significance  Level 
.10 
.05 
.01 


1006 

.0970 

.1003 

0486 

.0486 

.0498 

0103 

.0101 

.0102 

13 


3.   DESCRIPTION  OF  SIMULATION 

In  our  simulation  we  generated  deviates  from  three  families  of  distri- 
butions.  The  family  of  asymmetric  stable  distributions  has  been  previously 
used  by  Saniga  and  Miles  (1979)  to  investigate  the  power  of  several  goodness- 
of-fit  tests.   We  used  the  same  set  of  parameter  values,  a  =  1.0(.3)1.9  and 
8  ■  0(. 25) 1.00  they  used.   Mixed  normal  distributions  have  often  been  used 
for  such  tests,  and  we  have  included  a  family  which  is  a  composite  of  N(0,1) 
and  N(0,a),  a  ^  1,  and  another  set  which  is  a  composite  of  N(.5,l)  and 
N(0,a).   Pearson  distributions  considered  include  a  variety  of  shapes,  from 
'hear"  normal  to  U  and  J  shaped.   Discussion  of  the  procedures  used  to 
generate  pseudorandom  deviates  from  each  of  these  families  follows. 

The  sample  data  was  obtained  by  starting  with  one  or  more  uniformly  dis- 
tributed pseudorandom  deviates.   These  were  generated  using  the  IMSL  sub- 
routine GGUBS,  the  basis  of  which  is  discussed  in  Lewis,  Goodman,  and 
Miller  (1969). 

3.1  Generation  of  Asymmetric  Stable  Deviates  (Random  Stabilized  Standard 
Form) 

This  family  of  distributions  contains  as  a  special  case  the  normal 

distribution  (a  =  2 ,  8=0)  and  the  Cauchy  distribution  (a  ■  1,  8=0).   For 

a  <  2  the  distributions  have  infinite  variance,  which  makes  them  useful  for 

determining  the  ability  of  a  goodness-of-f it  test  to  detect  heavy  tailed 

distributions.   The  deviates  were  generated  using  the  program  given  in 

Chambers,  Mallows,  and  Stuck  (1976).   This  subroutine,  RSTAB,  used  one  deviate 

from  U[0,1]  and  one  exponentially  distributed  deviate  to  generate  one  RSSF 

deviate. 
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3.2  Generation  of  Mixed  Normal  Deviates 

Mixed  normals  of  the  form  (l-y)N(y..  ,o\)  +  y   N(u  .a.)  were  generated  using 

the  IMSL  routine  MDNRIS  to  convert  uniform  samples  into  standard  normal 

samples.   To  obtain  a  set  of  N  mixed  normal  variates  we  proceeded  as  follows: 

(i)  2N  uniform  random  variates  {u.}  were  generated  using  GGUBSr   (ii)  for 

each  i  =  1.....N,  MDNRIS  was  used  to  convert  u.  to  a  standard  normal  z.. 

1  1 

z.  was  then  transformed  to  a  normal  with  mean  u,  and  variance  a,,  or  with 
l  1  1 

mean  u„  and  variance  a„,  depending  on  whether  u.,„>Y,  or  u.,„  <  Y, 
2  2  x+N  -       i+N 

respectively. 

3.3  Generation  of  Pearson  Type  I,  II  Deviates 

The  generation  of  samples  from  Pearson  Type  I  and  II  distributions  was 

done  via  table  look  up  and  linear  interpolation  on  the  inverse  cumulative 

distribution  function.   Sufficient  entries  to  assure  four  significant  decimal 

places  in  the  final  answer  was  achieved  adaptively  using  numerical  integration. 

Before  discussion  of  the  precise  details  of  the  process,  we  digress  for  a 

discussion  of  the  Pearson  distributions,  and  particularly  types  I  and  II. 

Following  Johnson  and  Kotz  (1970),  the  Pearson  probability  density  function 

p(x)  is  given  by 

1_  dp_  _     a  +  x 

p  dx  2 

r        C   +  C.X  +  C-X 

o    1     2 

It  can  be  shown  that  a,  c  ,  x.,  and  c_  can  be  expressed  in  terms  of  non- 

o   1       I 

negative  parameters  B,  ,  B9,  and  the  variance  y„,  obtaining 

c  =  (462  -  331)(10B2  -  1231  -  18)_1u2 
a  =  c±   =/g^  (B2  +  3)(10B2  -  12B1  -  18)"1  /^  (7) 

c2  =  (2B2  -3B1  -6)(10B2  -  U^  -   18)"1  . 
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Type  I  distributions  are  characterized  by 

K:    =   4  Cl  (coC2)    <  °  ' 
while  Type  II  have  k  =  0  with  8,  =   0,   82  <  3.   For  Type  I  and  II  the  prob- 
ability density  function  is  of  the  form 

m        m 
p(x)  =  K(x  -  a^  x(a2  -  x)  c    .  (8) 

Generation  of  a  sequence  of  pseudorandom  Pearson  Type  I  and  II  deviates 
involves  the  following  steps:   (i)  generate  a  sequence  of  deviates   {u.}, 
from  U[0,1];   (ii)  transform  these  to  Pearson  deviates  by  finding  v.  so  that 

/    p(x)dx  =  u..   This  necessitates  being  able  to  efficiently  obtain  the 


ai 

inverse  CDF,  i.e.,  if  F(v)  =   /   p(x)dx,  then  we  need  v.  =  F   (u.).   We 


f 

[v)   =   /   p(x)dx, 


-1 
will  denote  F   by  G  in  order  to  simplify  notation. 

Representation  of  the  inverse  CDF,  G(s),  was  achieved  by  linear  inter- 
polation in  a  table  generated  by  numerical  integration  of  p(x).  The  adap- 
tive process  used  to  assure  four  decimal  place  accuracy,  i.e.,  magnitude  of 

-4 
the  error  less  than  . 5  x  10   is  described  now. 

The  error  in  linear  interpolation  between  points  (s.,G.)  and  (s .  . ,  G .  - ) 

1  2  i  "    i 

is  no  more  than  — (s .   -  s.)  (*  ,  where  G7  =   max     |G  (s)  |  . 

s.<s<s.  . . 
l-  -  l+l 

Since  the  intervals  will  be  small,  we  can  approximate 

G'(s.+1)  -  G'(s.) 

G„  by —   ,     and  then  using  the  fact  that 

2       s, x1  -  s. 

l+l    i 

G'(s)  =  (F_1(s))'  =  j^~   =  ^   ,  we  obtain 
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p(G.+1)    p(G.) 

s...  -  s. 
l+l    1 


Thus,  an  error  estimate  in  (s  ,  s .,,)  is  given  by 


p(G_)    p(G.) 


i+1 


s.,,  -  s. 
l+l    i 


(si+l  "  Si} 


p(G.+1)    p(G±) 


(s 


i+1 


Si> 


(9) 


Potential  problems  occur  if  p(G)  =  0,  as  may  happen  at  a  and  a  .   Since  we 

-4 
require  error  less  than  . 5  x  10   ,  it  is  clear  this  must  be  the  case  when 

I         i        -4 
G..,  -  G.  <.5  x  10   :  hence  if  p(G)  is  very  small  in  the  interval  (s. ,s. .,) 

we  have  an  alternative  scheme  for  accepting  an  interval,  one  which  came  into 

play  for  ranges  where  p  is  small. 

These  ideas  were  the  basis  for  adaptive  construction  of  a  suitable  table 

(s.,G.)  for  the  inverse  cumulative  distribution  function.   We  have 
l   l 

s  =  0,  G  =  a,.   We  describe  the  general  scheme  for  obtaining  (s .  , ,  ,  G..,) 

o        o     1  1+1    1+1 

given  (s.,G.)  and  note  how  an  estimate  of  G.M  is  generated  afterwards.   Given 
l   l  l+l 

an  estimate  for  G.,,,  the  value  of  As.  =  s . , ,  -  s.  is  obtained  by  numerical 
l+l  l    l+l    l 


numerical  integration  of 


J 

l 


Gi+i 


p(x)dx.   This  is  accomplished  by  the 


adaptive  quadrature  routine,  DCADRE,  from  the  IMSL  library.   If  m  or  m~  are 
negative,  subtracting  out  the  singularity  was  used  for  intervals  near  a.,  and 

a_,  respectively.   The  double  precision  version  of  DCADRE  is  used  and  an 

—6 
absolute  error  tolerance  of  10   is  requested.   The  routine  returns  As.  and 

-4 
the  error  estimate  (9)  is  computed.   If  it  is  less  than  . 5  x  10   ,  the 

result  is  accepted,  we  set  s.,1  =  s.  +  As.,  and  proceed  to  the  next  interval. 

-4  2 

Suppose  the  error  estimate,  E    >  . 5  x  10   .   Since  E    =  0((AG)  ),  and 

est  est 
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0(AG.)  =  0(As.),  we  obtain  a  new  estimate  for  AG.  by  taking  it  to  be 
11/  ^       1 

AGl//.25EXlO-4    ,    t 

est 

This  takes  the  error  based  on  the  new  AG.  to  approximately  midway  between  its 

-4 
current  value  and  .25  x  10   .   More  than  one  correction  of  this  sort  may  be 

required,  in  particular  initially  where  a  reasonable  estimate  of  AG  is  not 

o 

available. 

Once  an  interval  has  been  accepted  (or  rather,  a  point  (s.,G.))  we 
increment  the  interval  counter  i,  and  then  estimate  the  new  interval  size 


AG .  by  AGi-l  I .   /.475  x  10  4 

l    — - —   w +  1  J.   This  yields,  based  on  the  above  assumptions, 

V  »     est 

a  AG.  which  should  give  an  error  for  the  next  interval  which  is  midway  between 

-4 
the  previous  one  and  .475  x  10 

An  initial  value  for  AG  is  required.   Although  not  very  sophisticated, 

a2  "  al 

we  simply  take  AG  =   .  _._... —  ,  and  depend  on  the  adaptive  machinery  described 
o    1UUU 

above  to  decrease  it  to  meet  the  error  tolerance,  or  increase  successive 
intervals  as  required  for  efficient  representation. 

It  is  true  that  the  final  value  of  s.,  call  it  s  ,  should  be  equal  to  one. 
Because  of  the  numerical  integration,  this  is  never  exactly  the  case.   In 

the  worst  case,  3-,  *  .01,   3_  =  1.9,  we  obtained  the  final  value 

—8 
s  ~    .9999999845  an  error  of  about  1.55  x  10   .   In  order  to  avoid  any  prob- 
lems due  to  the  table  not  covering  [0,1]  exactly  the  simplest  procedure 
was  to  replace  each  computed  s.  by  s./s  ,  thus  distributing  the  error  over 

the  entire  interval  and  yielding  a  consistent  table.   Note  that  this  is 

-4 
well  within  the  error  tolerance  of  .5  x  10 

The  procedure  has  been  thoroughly  tested  for  its  efficiency  in  repre- 
senting the  inverse  CDF  as  well  as  for  accuracy.   For  the  particular  cases 
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in  which  we  were  interested,  the  inverse  CDF  was  represented  by  a  table  with 
no  more  than  729  entries.   This  occurred  for  the  case  8  =  .25,   B9  =  3.2, 
which  has  an  inverse  CDF  with  very  large  slopes.   For  more  gently  sloping 

inverse  CDF's  we  were  able  to  use  as  few  as  101  intervals,  as  in  the  case 

-4 
6  =  .01,   8  =  1.75.   Most  intervals  had  an  error  estimate  of  between  .4  x  10 

-4 
and  . 5  x  10   ,  which  shows  that  the  interval  sizing  process  we  have  used 

worked  quite  efficiently,  with  few  initial  estimates  being  rejected  for  being 

too  large,  without  also  resulting  in  intervals  much  too  small. 

The  routine  was  checked  against  the  published  tables  of  Johnson,  Nixon, 

Amos  and  Pearson  (1963)  for  many  values  of  the  parameters  and  at  most  of  the 

percentage  points.   With  a  few  exceptions,  where  a  difference  of  one  in  the 

fourth  decimal  place  was  noted,  the  results  check  exactly.   Generally  in 

these  cases  the  fifth  place  was  four  or  five  so  that  the  actual  error  was 

probably  well  within  our  tolerance. 
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4.   RESULTS 

The  results  of  the  simulation  are  summarized  in  Tables  5-12.   The  empi- 
rical power,  in  5000  replications,  of  the  Foutz  F  test  compared  with  that 

n 

for  the  Chi  squared  test  or  the  K-S  test  is  presented  as  a  percent  improve- 
ment in  the  probability  of  rejecting  the  null  hypothesis  that  the  distri- 
bution of  the  samples  is  the  standard  normal.   A  negative  entry  means  that 
the  power  of  the  F  test  was  smaller  than  that  for  the  Chi  squared  test  or 
the  K-S  test,  whichever  is  appropriate. 

The  simulation  has  revealed  that  the  F  test  is  better  than  the  Chi 

n 

squared  test  which  in  turn  is  better  than  the  K-S  test  when  the  true  distri- 
bution of  the  samples  is  heavy  tailed.   Many  such  distributions  are  included 
in  the  asymmetric  stable  family  as  well  as  the  family  of  mixtures  of  normals. 
For  the  mixed  normal  family,  if  the  two  normals  involved  in  the  mixture 

differ  in  the  means  the  K-S  test  performed  better  than  the  F  test  even  when 

n 

the  variances  differed.   We  now  discuss  in  more  detail  the  results  for  each 
of  the  three  families  of  distributions. 

4.1  Asymmetric  Stable  Family 

The  results  for  the  F  test  versus  the  Chi  squared  test  are  summarized 

n 

in  Table  5.   The  F  test  outperformed  the  Chi  squared  test  for  n  equal  to 

21  and  31.   When  n  ■  51,  as  a  +  g  increased  the  performance  of  the  F  test 

'  r  n 

deteriorated  as  can  be  seen  from  the  lower  right  part  of  Table  5.   Another 
general  observation  is  that  as  the  significance  level  is  decreased  the  improve- 
ment in  power  for  the  F  test  is  accentuated. 

n 

The  comparative  figures  for  the  F  test  versus  the  K-S  test  are  pre- 

n 

sented  in  Table  6.   Here  again  the  F  test  did  much  better  than  the  K-S  test; 

n 
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the  results  also  indicate  that  for  the  asymmetric  stable  family  the  Chi 
squared  test  has  a  higher  power  than  the  K-S  test. 

4.2  Mixtures  of  Normals 

The  mixed  normal  distributions  that  we  considered  were  of  two  basic 

types.   The  first  type  is  of  the  form  (l-y)N(0,l)  +  yN(0,a)  with  a  =  2,3,4 

and  y  =  -1>  «2,  .3,  1.0;  note  that  when  y   =  1.0  the  distribution  is  not  a 

true  mixture  but  N(0,o).   The  second  type  is  a  mixture  of  the  form 

(1-Y)N(.5,1)  +  yN(0,a)  with  a  =   3,    y  =    .2,    .3  and  a   =  4,  y   =  .2. 

The  F  test  did  significantly  better  than  the  Chi  squared  test  except 
n 

for  n  =  51  and  y  =  1-0  (see  Table  7);  in  the  latter  case  the  Chi  squared 

test  turned  out  to  be  the  better  of  the  two  tests. 

When  the  F  test  was  compared  with  the  K-S  test  the  F   test  turned  out 
n  n 

to  be  consistently  better  (Table  8). 

In  the  case  of  a  mixture  of  the  second  type, which  included  a  location 

shift  (Tables  9,  10), the  K-S  test  proved  to  be  superior  to  the  F  test 

n 

while  the  comparison  between  the  F  test  and  the  Chi  squared  test  appeared 

n 

to  be  inconclusive. 

4.3  Pearson  Family 

We  chose  ten  distributions  of  types  I  and  II  to  encompass  a  variety 

of  shapes  as  shown  in  figure  1;  the  standard  normal  is  superimposed  as  a 

dotted  curve  in  each  of  the  graphs.   The  comparison  of  the  F  test  versus 

n 

the  Chi  squared  test  proved  to  be  inconclusive.   However,  the  K-S  test 
appeared  to  have  a  higher  power  than  the  F   test  when  the  shape  of  the 
distribution  was  "near"  normal  such  as  the  ones  for 
(B1,B2)  =  (0,2.3),  (0,2.8)  and  (.25,3.2). 
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TABLE  5 


FOUTZ  F  TEST  VS  CHI  SQUARED  TEST 
n 

RANDOM  STABILIZED  STANDARD  DISTRIBUTIONS 


0.0 


0.25 


0.50 


0.75 


1.00 


n=21 


n=31 


n=51 


Sig 

Level 

a 

1.0 

i      102.0 

102.7 

101.6 

95.8 

96.9 

.10 

1.3 

96.1 

88.6 

69.1 

49.1 

11.7 

1.6 

72.0 

68.3 

59.5 

49.0 

29.2 

1.9 

54.3 

55.4 

55.4 

53.9 

55.8 

1.0 

120.2 

117.2 

121.6 

114.9 

110.3 

.05 

1.3 

97.9 

98.2 

71.2 

49.6 

-     1.4 

1.6 

71.2 

63.4 

58.6 

41.4 

20.1 

1.9 

49.8 

45.6 

51.1 

49.8 

44.2 

1.0 

256.0 

215.9 

238.0 

234.5 

233.8 

.01 

1.3 

174.8 

158.6 

131.8 

71.7 

-     8.0 

1.6 

118.5 

130.3 

107.9 

62.1 

21.7 

1.9 

81.6 

73.0 

99.0 

89.7 

78.7 

1.0 

37.9 

36.8 

38.1 

39.3 

37.0 

.10 

1.3 

35.9 

34.0 

28.4 

13.8 

-10.4 

1.6 

26.6 

26.2 

16.5 

11.1 

0.9 

1.9 

17.6 

19.5 

19.4 

13.4 

16.7 

1.0 

47.4 

44.6 

46.6 

50.2 

46.7 

.05 

1.3 

41.8 

38.9 

31.5 

11.7 

-18.6 

1.6 

28.8 

25.7 

13.0 

5.2 

-   7.6 

1.9 

12.0 

15.8 

13.3 

8.4 

14.4 

1.0 

74.7 

75.2 

73.8 

73.5 

70.2 

.01 

1.3 

55.8 

46.2 

31.2 

3.9 

-33.7 

1.6 

32.2 

30.7 

11.0 

-  6.0 

-21.3 

1.9 

14.5 

18.1 

10.1 

8.0 

1.0 

1.0 

11.4 

13.3 

11.1 

10.6 

11.4 

.10 

1.3 

9.6 

12.3 

6.1 

-  1.1 

-11.1 

1.6 

1.2 

-  0.5 

-  4.9 

-   7.3 

-13.9 

1.9 

-10.8 

-  8.9 

-11.0 

-11.1 

-11.2 

1.0 

17.4 

18.5 

17.1 

16.3 

17.1 

.05 

1.3 

12.6 

12.2 

5.0 

—  3.5 

-18.2 

1.6 

0.5 

-  3.5 

-  8.7 

-12.0 

-21.7 

1.9 

-14.3 

-12.9 

-15.1 

-16.3 

-15.9 

1.0 

30.5 

30.5 

31.2 

28.4 

25.8 

.01 

1.3 

14.1 

13.8 

3.1 

-12.2 

-35.2 

1.6 

-  5.4 

-  8.3 

-20.6 

-23.5 

-37.0 

1.9 

-24.2 

-24.2 

-24.1 

-24.2 

-26.  L 
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TABLE  6 

FOUTZ  F  TEST  VS  KOLMOGOROV-SMIRNOV  TEST  RANDOM 
RANDOM  STABILIZED  STANDARD  DISTRIBUTIONS 


Sig 


0.0 


0.25 


0.50 


0.75 


1.00 


LeveL 

a 

1.0 

120.2 

123.8 

110.0 

120.4 

115.2 

• 

.10 

1.3 

102.8 

100.6 

76.0 

43.5 

0.9 

1.6 

87.7 

77.2 

61.3 

49.6 

22.7 

1.9 

65.5 

61.7 

62.9 

65.1 

57.8 

1.0 

197.5 

196.0 

186.7 

185.9 

192.1 

.05 

1.3 

170.8 

151.2 

111.3 

68.1 

-     2.8 

n=21 

1.6 

131.3 

110.4 

89.9 

67.0 

30.4 

1.9 

83.7 

90.5 

87.5 

85.6 

74.4 

1.0 

601.5 

465.5 

530.4 

522.7 

526.8 

.01 

1.3 

396.4 

338.1 

245.0 

121.4 

-     0.6 

1.6 

305.7 

243.9 

165.2 

117.7 

45.6 

1.9 

140.1 

147.5 

176.2 

176.5 

163.8 

1.0 

82.1 

80.6 

83.3 

84.7 

81.6 

.10 

1.3 

83.6 

77.6 

56.7 

25.5 

-     9.2 

1.6 

63.8 

59.5 

42.6 

29.4 

7.3 

1.9 

49.5 

55.6 

48.2 

44.8 

40.3 

1.0 

150.9 

149.8 

145.6 

152.7 

151.5 

.05 

1.3 

147.3 

129.2 

88.1 

40.5 

-  13.0 

n=31 

1.6 

107.8 

98.9 

65.0 

41.1 

11.4 

1.9 

77.2 

86.7 

69.5 

72.5 

68.2 

1.0 

428.5 

429.5 

447.0 

399.0 

471.3 

.01 

1.3 

386.1 

316.4 

190.6 

76.4 

-   16.2 

1.6 

267.3 

210.3 

136.9 

79.9 

24.1 

1.9 

190.8 

186.2 

153.0 

146.8 

136.8 

1.0 

38.5 

42.7 

38.9 

37.5 

37.7 

.10 

1.3 

49.3 

49.4 

29.6 

8.4 

-  10.3 

1.6 

40.0 

35.7 

24.0 

11.4 

-     6.7 

1.9 

24.4 

31.4 

26.1 

22.2 

19.5 

1.0 

77.6 

81.2 

78.0 

78.3 

76.0 

.05 

1.3 

94.4 

85.2 

53.2 

14.7 

-   16.4 

1.6 

79.4 

61.4 

37.9 

17.0 

-     8.1 

n=51 

1.9 

49.8 

54.6 

50.6 

42.0 

35.4 

1.0 

285.4 

287.0 

275.2 

270.9 

262.6 

.01 

1.3 

287.6 

2  43.0 

138.7 

37.7 

-  26.1 

1.6 

196.2 

173.3 

91.2 

41.0 

-     7.9 

1.9 

130.4 

126.3 

115.4 

111.8 

99.6 
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TABLE  7 


FOUTZ  F  TEST  VS  CHI  SQUARED  TEST 
MIXED  NORMAL  DISTRIBUTIONS 


1.00 


0.30 


0.20 


0.10 


..     Sig 

-      N(0,a) 

Level 

- 

N(0,2) 

55.8 

69.3 

76.3 

40.7 

.10 

N(0,3) 

32.8 

86.8 

75.9 

44.7 

N(0,4) 

20.5 

101.4 

89.0 

75.7 

N(0,2) 

47.0 

46.4 

64.5 

23.7 

ii=21 

.05 

N(0,3) 

30.5 

71.3 

69.0 

34.7 

N(0,4) 

22.6 

91.1 

76.8 

48.1 

N(0,2) 

62.0 

262.1 

107.5 

88.6 

.01 

N(0,3) 

61.2 

109.8 

145.8 

86.3 

N(0,4) 

49.0 

118.3 

131.1 

125.0 

N(0,2) 

13.8 

46.9 

42.4 

22.0 

.10 

N(0,3) 

4.3 

56.7 

50.6 

51.5 

N(0,4) 

1.8 

55.8 

58.7 

39.2 

N(0,2) 

12.0 

38.3 

29.2 

28.6 

n=31 

.05 

N(0,3) 

3.1 

53.2 

48.4 

53.2 

N(0,4) 

1.0 

57.4 

61.5 

44.4 

N(0,2) 

1.2 

62.1 

46.0 

52.1 

.01 

N(0,3) 

0.6 

47.3 

103.1 

104.7 

N(0,4) 

0.1 

84.0 

95.5 

61.5 

N(0,2) 

-13.1 

23.6 

26.1 

9.7 

.10 

N(0,3) 

-  5.9 

19.3 

28.4 

25.7 

N(0,4) 

-   1.4 

26.2 

29.5 

35.6 

N(0,2) 

-16.2 

28.9 

45.5 

17.4 

n=51 

.05 

N(0,3) 

-  9.4 

29.2 

43.4 

38.0 

N(0,4) 

-   2.5 

34.1 

45.4 

39.7 

N(0,2) 

-26.9 

83.8 

54.8 

67.4 

.01 

N(0,3) 

-17.7 

41.6 

45.9 

67.9 

N(0,4) 

-   7.9 

24.5 

78.3 

51.4 
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TABLE  8 

FOUTZ  F  TEST  VS  KOLMOGOROV-SMIRNOV  TEST 
n 

MIXED  NORMAL  DISTRIBUTIONS 


Sig 


Y- 
N(0,a) 


1.00 


0.30 


0.20 


0.10 


N(0,2) 

60.5 

36.1 

31.8 

8.7 

.10 

N(0,3) 

52.6 

59.3 

42.9 

19.3 

N(0,4) 

38.8 

71.5 

54.0 

40.2 

N(0,2) 

87.4 

49.5 

36.4 

6.8 

n=21 

.05 

N(0,3) 

N(0,4) 

91.8 
74.4 

75.4 
90.2 

49.8 
66.6 

16.9 
43.7 

N(0,2) 

156.3 

94.4 

43.1 

15.8 

.01 

N(0,3) 

223.6 

112.3 

71.0 

41.8 

N(0,4) 

198.1 

107.0 

101.2 

26.6 

N(0,2) 

42.6 

19.3 

27.8 

13.9 

.10 

N(0,3) 

31.3 

52.9 

38.8 

24.1 

N(0,4) 

18.5 

66.2 

52.9 

21.1 

N(0,2) 

69.3 

29.9 

27.8 

20.6 

n=31 

.05 

N(0,3) 
N(0,4) 

60.8 
38.7 

65.7 
101.4 

52.9 
61.9 

36.8 
36.2 

N(0,2) 

135.8 

52.9 

67.3 

55.3 

.01 

N(0,3) 

183.9 

107.5 

120.3 

66.0 

N(0,4) 

129.8 

163.9 

152.2 

72.1 

N(0,2) 

22.4 

26.8 

23.6 

9.7 

.10 

N(0,3) 

9.6 

47.8 

39.4 

21.2 

N(0,4) 

3.1 

66.4 

51.1 

36.3 

N(0,2) 

45.7 

30.0 

39.9 

9.0 

n=51 

.05 

N(0,3) 
N(0,4) 

24.5 
10.6 

74.0 
110.7 

52.2 
73.3 

29.9 
52.3 

N(0,2) 

117.6 

61.9 

81.1 

20.3 

.01 

N(0,3) 

96.0 

215.7 

83.3 

39.1 

N(0,4) 

53.0 

170.9 

164.5 

86.7 

1 
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TABLE  9 


FOUTZ  F  TEST  VS  CHI  SQUARED  TEST 
n 

MIXED  NORMAL  DISTRIBUTIONS 


MIXED  NORMAL 


SIGNIFICANCE  LEVEL 
10        .05 


,01 


n=21 

0 .  7GxN  (0 . 5 ,  l)+0 .  30xN  (0 .0,3) 
0.80xN(0.5,l)+0.20xN(0.0,3) 
0.80xN(0.5,l)+0.20xN(0.0,4) 

24.0 
14.4 
17.2 

-  0.7 

-  9.2 

-  3.2 

-20.8 
-27.2 
-21.1 

n=31 

0.70xN(0.5,l)+0.30xN(0.0,3) 
0.80xN(0.5,l)+0.20xN(0.0,3) 
0.80xN(0.5,l)+0.20xN(0.0,4) 

9.6 
3.3 
3.7 

-   6.8 
-12.4 
-10.7 

-28.6 
-33.3 
-30.2 

n=51 

0.70xN(0.5,l)+0.30xN(0.0,3) 
0.80xN(0.5,l)+0.20xN(0.0,3) 
0.80xN(0.5,l)+0.20xN(0.0,4) 

25.1 
-   2.3 

-24.1 
-23.6 
-26.3 

-43.6 
-48.6 
-45.5 
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TABLE  10 


FOUTZ  F  TEST  VS  KOLMOGOROV-SMIRNOV  TEST 
n 
MIXED  NORMAL  DISTRIBUTIONS 


MIXED  NORMAL 


SIGNIFICANCE  LEVEL 
.10        .05        .01 


n=21 

0.70xN(0.5,l)+0.30xN(0.0,3) 
0 . 80xN (0 . 5 , l)+0 . 20xN (0.0,3) 
0.80xN(0.5,l)+0.20xN(0.0,4) 

-17.0 
-27.3 
-25.5 

-28.9 
-38.2 
-32.5 

-33.7 
-40.7 
-35.1 

n=31 

0.70xN(0.5,l)+0.30xN(0.0,3) 
0.80xN(0.5,l)+0.20xN(0.0,3) 
0.80xN(0.5,l)+0.20xN(0.0,4) 

-25.4 
-34.8 
-34.6 

-34.5 
-45.5 
-40.2 

-43.4 
-50.7 
-45.9 

n=51 

0.70xN(0.5,l)+0.30xN(0.0,3) 
0.80xN(0.5,l)+0.20xN(0.0,3) 
0.80xN(0.5,l)+0.20xN(0.0,4) 

-28.7 
-48.7 
-43.3 

-47.6 
-57.4 
-51.3 

-56.7 
-64.9 
-61.0 
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TABLE  11 


FOUTZ  F  TEST  VS  CHI  SQUARED  TEST 
PEARSON  DISTRIBUTIONS 


i 

SIGNIFICANCE 

LEVEL 

si 

82 

.10 

.05 

.01 

0.0 

2.30 

3.4 

-10.5 

-20.6 

0.0 

2.80 

21.1 

-  0.9 

-12.4 

0.01 

1.75 

12.8 

41.1 

57.8 

0.01 

1.90 

-  9.1 

-  2.4 

18.9 

n=21 

0.25 
0.50 

3.20 
3.00 

16.6 
17.5 

-  0.4 
5.8 

-  2.3 

-  9.3 

1.00 

3.40 

-12.5 

-16.2 

-15.8 

1.00 

3.60 

9.1 

-31.4 

-30.9 

1.00 

3.80 

16.1 

-17.6 

-24.8 

2.00 

5.50 

-16.1 

-14.9 

-20.7 

0.0 

2.30 

-  6.6 

-19.0 

-22.5 

0.0 

2.80 

10.7 

-12.1 

-  5.4 

0.01 

1.75 

-  5.3 

39.4 

83.5 

0.01 

1.90 

-18.7 

-  6.9 

26.9 

a-31 

0.25 
0.50 

3.20 
3.00 

0.7 
-  1.1 

-11.8 
-  0.6 

-  3.6 
-12.2 

1.00 

3.40 

-25.5 

-23.5 

-23.4 

1.00 

3.60 

-  5.7 

-41.4 

-41.7 

1.00 

3.80 

0.0 

-26.9 

-32.9 

2.00 

5.50 

-29.5 

-19.4 

-26.5 

0.0 

2.30 

30.4 

-39.0 

-19.7 

0.0 

2.80 

34.9 

-40.0 

-17.3 

0.01 

1.75 

6.6 

68.1 

132.9 

0.01 

1.90 

-11.8 

-15.8 

20.7 

n-51 

0.25 
0.50 

3.20 
3.00 

33.3 
2.5 

-10.8 
-19.9 

-15.6 
-20.2 

1.00 

3.40 

-37.4 

-35.7 

-35.6 

1.00 

3.60 

-13.6 

-53.7 

-56.7 

1.00 

3.80 

-  4.7 

-40.6 

-46.5 

2.00 

5.50 

-38.5 

-28.7 

-38.0 
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TABLE  12 


FOUTZ  F  TEST  VS  KOLMOGOROV-SMIRNOV  TEST 
n 

PEARSON  DISTRIBUTIONS 


.10 


SIGNIFICANCE  LEVEL 
.05 


,01 


0.0 

2.30 

-22.8 

-25.0 

-26.2 

0.0 

2.80 

5.4 

-18.8 

-   7.2 

0.01 

1.75 

-  9.4 

20.6 

66.0 

0.01 

1.90 

-33.3 

-15.3 

16.5 

n-21 

0.25 
0.50 

3.20 
3.00 

-  8.1 
9.9 

-14.8 
21.7 

-  4.1 
34.9 

1.00 

3.40 

84.8 

94.4 

95.9 

1.00 

3.60 

52.9 

64.5 

76.7 

1.00 

3.80 

38.4 

44.0 

59.4 

2.00 

5.50 

62.4 

74.7 

70.4 

0.0 

2.30 

-28.9 

-36.8 

-36.9 

0.0 

2.80 

11.2 

-20.8 

-13.0 

0.01 

1.75 

-14.1 

27.6 

89.1 

0.01 

1.90 

-36.7 

-15.4 

19.4 

n-31 

0.25 
0.50 

3.20 
3.00 

-12.9 
-  0.6 

-14.0 
25.0 

-  9.3 
31.3 

1.00 

3.40 

100.0 

134.3 

131.7 

1.00 

3.60 

55.5 

79.4 

91.8 

1.00 

3.80 

39.9 

46.2 

65.3 

2.00 

5.50 

63.3 

91.0 

100.0 

0.0 

2.30 

-  9.1 

-48.6 

-40.0 

0.0 

2.80 

-  4.9 

-38.9 

-35.8 

0.01 

1.75 

-18.4 

62.0 

171.8 

0.01 

1.90 

-37.4 

-22.9 

26.6 

n=51 

0.25 
0.50 

3.20 
3.00 

-   8.1 
0.0 

-  9.6 
25.9 

-13.8 
47.0 

1.00 

3.40 

169.9 

225.8 

245.5 

1.00 

3.60 

64.9 

138.9 

154.4 

1.00 

3.80 

42.1 

69.0 

85.2 

2.00 

5.50 

103.0 

151.0 

164.7 
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5.   CONCLUDING  REMARKS 

The  superior  performance  of  the  Foutz  F  test  in  detecting  certain  types 
of  deviations  from  the  hypothesized  distribution  leads  to  several  more  prob- 
lems to  be  considered.   Of  primary  importance  is  the  generation  of  percentage 

points  for  the  distribution  of  F  for  various  values  of  n.   The  intractability 

n  J 

of  the  problem  of  obtaining  the  exact  distribution  requires  an  empirical 
approach  to  finding  a  correction  to  the  asymptotic  approximation  given  by 
Foutz. 

Since  the  test  is  also  applicable  to  p-variate  distributions,  an  investi- 
gation of  ways  to  obtain  the  statistically  equivalent  blocks,  and  then  the 
probability  content  of  them,  at  least  for  p  =  2,  is  to  be  considered.   The 
problems  of  obtaining  these  blocks  and  their  contents  becomes  increasingly 
complicated  in  higher  dimensions. 

The  authors  wish  to  acknowledge  extensive  use  of  the  IBM  360/67  computer 
system  at  the  Naval  Postgraduate  School,  on  which  all  computations  were 
performed. 
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APPENDIX  1 

C 

C      THIS  PROGRAM  GENERATES  DEVIATES  ?N  RANDOM  STABILIZED  STANDARD  FORM 

C      BLCCKS  GF  23,  30  ♦  ANC  50  DEVIATES  ARE  CONSIDERED  BY  GENERATING  53, 

C      THEN  CONSIDERING  THE  PIRST  20,  THE  NEXT  30,  AND  ALL  50. 

C      THE  HYPOTHESIS  THAT  THEY  CAME  FROM  A  NORMAL  DISTRIBUTION  WITH 

C      MEAN  AND  VARIANCE  AS  GIVEN  IN  INPUT  IS  TESTED. 

C 

C  THE    TESTS    PERFORMEC    ARE   THE    CHI    SQUARED,    KOLMOGOROV-SMIRNOV ,    AMD 

C  FOUTZ    TEST    AT    THE    CCNFIOENCE    LEVELS    OF    10,    5,     AND    1    PERCENT. 

C  THE    INPUT    VALUES    ARE    MEAN,     VARI ANCE, I  SEED ,    AND    THE    NUMBER    OF    REPLI 

C  CATIONS. 

C 

IMPLICIT    REAL*8    (A-Z) 

INTEGER*4   T, IER, N , I ,NP1,N1 , N2 ,N3, 1 SEED,NR ,NSTOP, IDF , IQ  ,ISEED1 

1  ,NCEL(3I,NST(3) ,N SN P (3 ) ,NCHI (3 ,3 ) ,NCHIN (3,3 ) ,  NKS (3,  3 ) , 

2  NKSN(2,3),NFNN(2,2),NFN(3,3) ,NSR 
DIMENSICf^    RD(551  .CHSQUI2)  , 

1    K0LSMR(2)  ,<=0UTZL(2)  ,FNTST(3,3)  ,  ZALF  (3  )  ,  RN(5  5  ) 
REAL*4    CELLS(50),CCMP(50),Q,CHI,R<  55) ,POIF( 6) , QTST(  3) ,RNN(55) , 

1    STST<3,3) 
DATA    ST  ST/.  26473,.  2 1 756,.  16<5  5  S,.  2940  8,.  241 70,.  18  841,.  35241, 

1    .28987, .22604/ 
DATA    NCEL/8,10,10/,NST/l,21,l/,NSMP/20,30,50/,QTST/.l,.05, .01/ 

1  ,CHISQU/»CHI    SQUAS'RED    TE  ST  • /,KOLSMR/»  KOLMOG-S*  ,•  MI  R    TEST*/ 

2  ,FOUTZL/«  FOUTZ    TE •,' ST» /, ZALC1, ZALC2/ .24306900, . 36787900/ 

DATA    WT1,MEAN1.VAR1,WT2,MEAN2,VAR2/.8D0,0.D0,1 .DO, .2 DO, 0 .DO ,3 .00/ 

3  , ZALF/ 1.28155 CO, 1.6448 500,2.3 263 500/ 
C 

C  THE    NEXT    4  CARDS    PUT    IN   THE    EMPIRICAL    CRITICAL    VALUES    FOR   THE 

C  FOUTZ    DISTRIBUTICN. 


C 


DIMENSION    FNEWP<3,3) 

DATA    FNEMP/. 4271400,. 41903DO,. 4)8160  0, 

1  .44865C0, .4355300, .4211600, .48659D0, 

2  .46579D0,.44487D0/ 

REAL*4  ALPHA, BETA, BPRIME 
COMMON/NCRPRM/MEAN,VARI,SVARI 
EXTERNAL  NCRM,UNIF 
T  =  1 

PIB2  =  CATAN(1.DC)*2.D0 
100  PRINT  1 

READ  2,MEAN,VARI 

IF (MEAN +VARI.LE.C.C0)  STOP 

2  F0RMAT(2E5.0> 
SVARI    =    DSORT(VARI) 
WRITE<4,14) 

14    FORMAT(lHl) 
1    FORMAT ( «+INPUT    MEAN    AND    VARIANCC,    FORMATC 2S5.0I • 1 

PRINT    17 

READ    12, ALPHA, BETA 
17    FORMAT!  •■*•    INPUT    ALPHA    AND    BETA,    FORMAT(  2F5. 0)  •  ) 
12    FORMATC5F5.0) 

PF  INT    3 

READ   5,ISEED,NP 

ISEED1    =    ISEEC 

BPRIME    =    BETA 

IF(ALPHA.EC.O. )GC    TC    110 

PHIZ    =    -PIB2*BETA*(1.    -    ABS(1.    -    ALPHA )) /ALPHA 

BPRIME    =    DTAN(PIB2*<1.    -    ALPHA  )  )*OTAN( PHI Z*ALPHA) 
110    CONTINUE 

DSEED    =    DFLOAT(ISEEC) 
5    FORMAT(  110,15) 

3  FORMAT( '+INPUT    SEED    AND    NUMBER    QF    REPLICATIONS,    F0RMAT( 110 , 1 5 ) • ) 
CO    205    10*1,3 

DO    205    M  =  l,3 
NCHKNl.IQI    =    C 
NCHIN(Nl.IQ)    =    0 
NKS(N1,IC)    =    0 
NKSN(N1,I0)    =    C 
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NFNCN1, IC I    =    0 

FNTSTCNl.ICI     =    ZAIC1*ZALF(IC)/DSCRT(DFL0AT(NSMP<N1I+1  ))    ♦    ZALC2 
C 

C  THE    FOLLCWING    CARC    PUTS    IN    THE    EMPIRICAL    VALUES    FOR    THE 

C  FOUTZ    TEST 


C 


FNTST(N1,IQ)    =    FNEMP<N1,IQ! 

205    NFNN(Nl.IO)    =    0 

DO   300    1=1, NR 

CALL    RANSRVC 50  ,RN  ,RD  ,OSSED, ALPHA ,BPRIME I 

CO    250    Nl=l,3 

N2    =    NST(Nl) 

N3    =    NSMP(Nl) 

CALL    VSRTAD(RD<N2),N3» 

CALL    VSR7AD(RN(N2) ,N3) 

DO    210    NSR=1,N2 

RNN(NSR)=    RN(NSR    +    N2-1) 
210    R(NSR)     =    RD(NSR+N2-1» 

IDF    =    0 

CALL    GFIT(UNIF,NCEL<N1),PNN,N3,CFLLS,C0MP,CHI  ,  IDF,Q,IER) 

CO    241    IC=1,3 

IF(0,LT.CTST(IOnNCHI(Nl,IQ)     =NCHI(N1,IQ)    +    1 

241  CONTINUE 
CALL    NKS1(UNIF,RNN,N3,PDIF,I5R) 
CO    242    IC    =    1,3 
IF(PDIF<1)  ,GT.STST(M,I0))NKS(N1,IC)    =    NKS(Nl.IQ)    +    1 

242  CONTINUE 
CALL    FCUTZ(UNIF,PNN,N3,Q) 
DO    243    IC-1,3 
IF(Q.GT.FNTST(N1,IQ) )NFN(N1,IQ>    =    NFN(N1,I0)    ♦    1 

243  CONTINUE 
IDF    =    0 

CALL    GFIT(N0RM,NCEL<N1),R,N3, CELLS, COMP,CHI ,IDF,Q,IEP) 
00    244    IC    =    1,3 
IFC0.LT.0TST(IQMNCHIN(N1,I0)     =    NCHIN(N1»I0)    +    1 

244  CONTINUE 
CALL    NKSl(N0RM,R^3,P0IFtIER) 
DO    245    10=1,3 
IF(PDIF(l).GT.STST(M,IQnNKSN(Nl,IQ)    =    NKSN<N1,IQ)    ♦    1 

245  CONTINUE 
CALL    F0UTZ(NCRN,R,N2,Q) 
DO    246    IC=1,3 
IF(Q.GT.FNTST(M  ,ICMNFNN<N1,IC)    =    NFNN(N1,I0)    ♦    1 

246  CONTINUE 
250  CONTINUE 
300    CONTINUE 

ISEED    =    CSEED 

WRITS<4,7»MEAN,VAFI,NP,!SFED1,ISEEC 
WRITE<4,11)ALPHA,BSTA 
11    FORMAT(//»    RANDOM    STABLE    STANDARDIZED    FORM    (RSSF),    ALPHA    =»,F5.2, 
1    •  ,    BETA    =«,F5.2//> 
WRITE(4,€)CHISQU,((NCHI(N1,IQ) , 

1  Nl=l,3)  ,IQ=1,3>  ,UNCHIN(M,IC>, 

2  N1=1,3),I0=1,3> 
WRITE(4,8)K0LSMR,<<NKS<N1,IQ),N1=1,3I,IQ=1,3) , 

1    UNKSNCM  ,IC)  ,N1=1,3),IQ=1,3) 

WRITE(4,8)F0UTZL,(<NFN(N1,IQ) ,N1  =  1 ,3 ) ,10=1 ,3 ) , i ( NFNN< M ,1  Q) ,N1=1, 
1     ),    10=1,3) 

7  FORMATC//'    THE    VALLES    OF    MEAN     (U)     AND   VARIANCE     (VI    ARE  ■ , 2F12.5/// 
9    ■    THE    RESULTS    OF    «,I5, •    P5PL ICATI CNS,    STARTING    WITH    SEED', 

1  I12//33X, 'ENDING    kITH   NEXT    SEECNI12//) 

8  FORMAT(///20X,2A8//20X,«20    PTS', 

2  6X,«30    PTS'.eX.'SO    PTS  •// 10X,  •  10*  •, 31 12/»    CONTROL •, 2X, •    5£«, 

3  3I12/10X,'    l*«,2I12///«       RSSP       • ,  2X  ,  •  10?1  ,3!  12/  •         VS',5X,f    5*», 

4  3112/*    N    <U,V)«,2X,»    1SS3I12) 

WRITE  (4,  9)  <<FNTST<M,IQ),N1  =  1,3),IC=1,3) 

9  FORMAT(/«0  7FIS    WAS    BASEC    ON    THE    FOLLLOWING   TEST    VALUES1, 
1    /(13X,3F12.5)  ) 

GO   to    100 
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END 


SUBROUTINE    RANSRV<ND,RAN,R, DSESD, ALPHA ,BPRIME) 
IMPLICIT    REAL*8     (A-y,0-Z) 
DIMENSION    P(l)  ,RAM1) 
REAL*4    RN( 100), ALPHA, BPRI ME, R STAB, W 
CALL   GGUBS(DSE6C,2*NC,RN) 
DO    200    NMltNO 
RAN(NN)     =    RN(NN) 
W    »   -BtOGUVOO    -    PN(NN+ND)) 
R(NN)    -    RSTAB( ALPHA, BPRI ME , RN ( NN) ,W) 
200    CONTINUE 
RETURN 
EHD 


SUBROUTINE    NCRMXtPJ 
REAL*8    MlANtVAPltSVARl 

COMMON/NCR  PR  l*/fEAh ,VARI ,SVARI 

T    =    (X    -   MEANI/SVARI 

P    =    .5*£PFC(-T*. 7071068) 

RETURN 
END 

SUBROUTINE   UNIF(X,P) 

P    =    X 

RETURN 

END 

FUNCTION  RSTAB(ALPHA,8PRIME,U,W) 
C      RSTA3  IS  A  RANDON  STABLE  STANDARDIZED  FORM  (WHATEVER) 
C 

C      ARGUMENTS 

C  ALPHA       =         CHARACTERISTIC    EXPCNENT 

C  BPPIME    =      SKEGNESS    IN    REVISED    PARAMETERIZATION 

C  U  UMFCPf    VARIATE    ON    (0,1) 

C  W  EXPCNENTIALLY    DISTRIBUTED   VARIATE 

DOUBLE    PRECISION    CA.CB 

DATA    P I BY2/1. 57079633/, PI BY4/. 785398 163/, THR1/. 99/ 

EPS    «    1.    -    ALPHA 
C  COMPUTE    SOME    TANGENTS 

PHIBY2     =    PIBY2*(L    -    .5) 

A    =    PHIBV2*TAN2(PHI8Y2) 

BB    =   TAN2(EPS*PHIEY2) 

B    =   EPS*FHIBY2*BB 

IF(EPS.GT.-.99)TAU    =    BPRIME / ( TAN2 (EPS*PIBY2) *P IBY2) 

IF(EPS.LE.-.99)TAL    =    BPRIME*PIBY2*EPS*M.    -E°S  )*TAN2(  ( 1.    -    EPS)* 
1    PIBY2) 
C  COMPUTE    SOME    NECESSARY    SUBEXPRESSIONS 

C  IF    PHI    NEAR    PI    BY    2,    USE    DOUBLE    PRECISION. 

IF(A.GT.THR1)G0    TO    50 
C  SINGLE    PRECISION 

A2    =   A**2 

A2P    =    1.    +   A2 

A2    =    1.    -    A2 

B2    =   B**2 

B2P    =    1.    +    B2 

B2    =    1.    -    82 

GO    TO    IOC 
C  COUBLE    PRECISION 

50    DA    «    CBLE(A)**2 

DE    -    DBLE(B)**2 

A2    =    l.CO    -    DA 

A2P    -    l.DO    +    DA 

B2   *    1.    -   OB 

C2P    =    1.    +   DB 
C  COMPUTE    COEFFICIENT 
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100    Z    -    A2P*(B2    ♦    2.*PHIBY2*3B*TAU)/(W*A2*82P) 
COMPUTE    THE    EXPONENTIAL-TYPE    EXPRESSION 
ALOGZ    =    ALOG(Z) 

D    -    D2(EPS*ALCGZ/(1.    -   EPS ))*( ALOGZ/ ( 1.    -   EPS)) 
COMPUTE    STABLE 

RSTAB    =    (1.    +    EPS*D)*2.*((A    -    B)*(l.    +    A*B)    -    PHI  BY2*TAU*BB*(  B*A2 
1    -    2.*A)  )/<A2*E2P)    ♦   TAU*n 
RETURN 
END 

FUNCTION    D2(Z) 
EVALUATE    (EXP(X)    -    1)/X 
DOUBLE    PPECISICN    Fl  ,  F2 ,01 , Q2 , C3, PV ,ZZ 

DATA    P1,P2,Q1,Q2,03/.8400668  5  2  5364£3  23903,.2')00111415  89964569D2, 
1    .168013370507926648C4,.1800133704C7390023D3, l.DQ/ 
THE    APPROXIMATION    1801    FOP    HART    ET    AL    (1968, P213) 
IF(ABS(Z).GT.0.1 )GC    TO    100 
ZZ    =    Z*Z 
PV    *    pi    +    ZZ*P2 

C2    =    2.C0*PV/(Q1    4    ZZ*(Q2    ♦    ZZ*Q3)    -    Z*PV) 
RETURN 
100    D2    =    (EXP(Z)    -    1.1/Z 
RETURN 
EMO 


FUNCTION    TAN(XARG) 

TANGENT    FUNCTION 

LOGICAL    NEG, INV 

DATA    P0*P1,P2, 00,01,02/    . 129221035E3 ,-.887662377El, .5286444565-1, 
1.1645  29332E3,-.45l!20561E2, 1./ 

THE    APPRCXIMATICN    4283    FROM    FAPT    ET    AL     (1968,     P.    251) 

DATA    PIBY4   / .7852*8163/ ,PI BY2/1. 57C79633/ , PI /3 .14159265/ 

NEG    =    .FALSE. 

INV    =    .FALSE. 

X    =    XARG 

NEG    =    X.LT.O. 

X    =    ABS(X) 

PERFORM    RANGE    REDUCTION    IF    NECESSARY 

IF(X.LE.FIBY4)GC    TC    50 

X    =    AMOD(X,PI) 

IF(X.LE.FIBY2)G0    TO    30 

NEG    =    . NCT.NEG 

X    ■   PI    -    X 
30    IF(X.LE.FIBY4)    GC    TC    50 

INV    =    .TRUE. 

X    =    PIBY2    -    X 
50    X    =    X/PIEY4 

CONVERT    TO    RANGE    CF    RATIONAL 

Y  v     —     y  &  v 

TAN    =    X*(PO    +    XX*(F1    ♦    XX*P2)»/(C0    ♦    XX*(Q1    +    XX*Q2l) 

IF(NEG)TAN    =   -TAN 

IFdNVJTAN    =    l./TAN 

RETURN 

END 


FUNCTION   TAN2(XARG) 
C  COMPUTE    TAN(X)/X 

C  FUNCTION    DEFINED    CNLY    FCR    ABS ( XARG) . LE. PI    BY    4 

C  FOR    OTHER    ARGUMENT    RETURNS    TAN(X)/X,    COMPUTED    DIRECTLY 

DATA    P 0,P1,P2, CO,  CI,  C2/. 12 9221035E3,-. 88766237 7E1, .5286444563-1, 
1     .16452<3332E2,-.4fl220561S2,l./ 
C  THE    APPRCXIMATICN    4283    FROM    HAPT    ST    AL     (1968,     P.    251) 

DATA    P I BV4,PIBY2,PI/. 785398 16 3, 1.57079633,3.14 15926 5/ 

X    *    ABS(XARG) 

IF(X.GT.PIBY4)GC    TC    230 

X    =    X/PIBY4 
C  CONVERT    TO    RANGE    CF    RATIONAL 
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c 

C       THIS  PRCGRAM  GENERATES  MIXED  NORMAL  DEVIATES   WT1*N ( MEAN1 , V4P 1 )  ♦ 

C       (1  -  WTl)*N(MEA|v2,VAR2)  . 

C  BLOCKS    OF    20,     30,    AND    50    DFEVIA^ES    ARE    CONSIDERED    BY    GENERATING 

C  50,    THEN    CCNSICEFINC    THE    FIRST    20,    THE    NEXT    30,    AND    ALL    50. 

C  THE    HYPOTHESIS    TMT    THEY    CAME    FROM    A    NORMAL    DISTRIBUTION    WITH 

C  MEAN    AND    VARIANCE    AS    GIVEN    IN    INPUT    IS    TESTED, 

C 

C  THE    TESTS    PERFORMED    ARE   THE    CHI    SQUARED,    KOLMOGOROV-SMIRNOV,    AND 

C  FOUTZ    TEST    AT    THE    CCNFIOFNCF    LEVELS    OF    10,    5,    AND    1    PERCENT, 

C  THE    INPUT    VALUES    APE    MEAN,    VAR IANCE, ISEED,    AND    THE    NUMBER    OF    REPLI 

C  CATIONS. 

C 

IMPLICIT    REAL*8    (A-ZI 

INTEGER* A    T,IER,N*I  ,NP1 ,N1 ,N2  ,N3 ,1  SEED, NR.NSTOP, I  OF, IQ,ISEE01 

1  ,NCEL(3),NST(3),NSr'P(3),NCHI<3,3),NCHlN(3,3),NKS(3,3)  , 

2  NKSN(3,3) ,NFNN(3,3) ,NFN<3 , 3 ) , NSR 
DIMENSION    RD(55> ,CHISQU(2) , 

1    KCLSMR(2I,FGUTZL<2),FNTST(3,3),ZALF(3),RN<55) 
REAL*4    CELLS (50 1 ,CCMF(50) , Q, CHI  ,R < 55  ) , PDIF<6 ) , OTST (3 ) , RNN ( 5  5) , 

1    STST(3,=) 
CATA    STST/. 26473,. 21 756,. 1695 9, ,29408,. 24170,. 18841,. 35241, 

1    . 28987, .226C4/ 
CATA    NCEL/8,10,10/,NST/l,21,l/,NSMP/20,30,50/,QTST/.l,.05,.01/ 

1  ,CHISQL/»CHI    SCUA* ,«RE0   TEST1/, KOLSMR/ 'KOLMOG-S • , 'MIR    TEST'/ 

2  ,FOUTZL/'FOUTZ    TE •,• ST' /,ZALC1,  ZALC2/ .24306900, .36787900/ 

DATA    WT1,MEAN1,VAR1,WT2,MEAN2,VAP2/.8D0,0.D0,1.D0,.2D0,0,D0,3.D0/ 

3  ,  ZALF/ 1.281 55  CO, 1.6 448 500,2.3  2635 DO/ 
C 

C  THE    NEXT    4   CARDS    PUT    IN    THE    EMPIRICAL    CRITICAL    VALUES    FOR   THE 

C  FOUTZ    TEST. 


C 


i! 


ClfSNSION    FNEMP(3,3) 

DATA  FNEMP/. 4271400,. 4 190 3D0,. 40 816D0, 

1  .4486500, .4355300,  .4211600, .486 59D0 , 

2  ,A657CD0,,44487D0/ 

:OMMQM/NCRPRM/^EAh ,VARI ,SVARI 
JXTERNAL    NORM,UNIF 
T    =    1 
100    PRINT    1 

READ  2,MSAN,VARI 
IF(MEAN+VARI,LE.O.CO)STOP 

2  FORMAT(2E5.0) 
SVARI    =    CSCRT(VAPI) 
WRITE  (4, 14 1 

14    FORMAT (1H1) 
1    FORMAT ('+INPUT    MEAN    AND    VARIANCE,     FORMAT  (  2E5.0  )  •  ) 

PRINT    17 

READ    12,WT1,MEAN1,VAR1,MEAN2,VAP2 

WT2    s    1    CO    -   WT1 
17    FORMATS*    INPUT    V»T1  ,MEAN1 ,  VAR1  ,MEAN2  ,VAR2  ,    FORMAT  C5F5  .0  )•  I 
12    FORMAT  ( 5  F5.0) 

PRINT  3 

READ  5, ISEED, NR 

ISEED1  =  ISEED 

DSEED  =  CFLOAT(ISEEC) 
5  FORMAT (110,15) 

3  FORMAT( '+INPUT  SEEC  AND  NUMBER  OF  REPLICATIONS,  FORMAT( 110 ,1 5)  •  ) 
CO  205  10*1,3 

DO  205  Nl=l,3 
NCHI(N1,IC)  =  0 
NCHIN(N1,IQ)  =  0 
NKS(N1, IC)  -    0 
NKSN(N1,I0)  =  0 
NFN(N1,IQ)  =  0 

FNTST(N1,IC)    =    ZALC1*ZALF(I0)/DS0RT(DFL0AT(NSMP(N1)+1) )    ♦    ZALC2 
FNTST(N1,IC>     -    FNENFCMtlQ) 
205    NFNN(N1,IQ)    =    0 
DO    300    I*1,NR 
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CALL    MIXNRM(50,RN,RD,DS5FD,WT1,MEAN1 , VAR1 , ME AN2 , VAR2 ) 
CO    250    M  =  l,3 
N2    =   NSTCNU 
N3    =    NSMP(Nl) 
CALL    VSRTAD(RD<N2),N3) 
CALL    VSRTAD(RN(N2),N3) 
CO    210    NSR=1,N3 
RNN(NSR)*    RN(NSR    4    N2-1) 
210    R(NSR)     =   RD(NSR+N2-1) 
IDF    =    0 

CALL    GFIKUNIF  ,NCEL<  Nl  I  ,RNN,  N3  ,CELLS  ,  COMP,CHI ,  IDF,Q, IERI 
CO    241    10=1,3 
IFCO.LT.CTSTUC)  )NCH(N1,!Q)  .=   NCHKN1,I0)    +    1 

241  CONTINUE 

CALL    NKSHUNIF,RNN,N3,PDIF,T  =  R) 

DO    242     IC    =    1,3 

IFfPDIFU)  .GT.STST<M,IQ))NKS(M,IC)     =    NKS(N1,IG)    +    1 

242  CONTINUE 

CALL  FCLTZ(UNIF,PNN,N3,Q) 

DO  243  10*1,3 

IF(Q.GT.FNTST<M,IC)  )NFN(N1,  IQ)  =  NFN(N1,IQ)  ♦  1 

243  CONTINUE 
IDF  =  0 

CALL  GF IT( NORM t NCE L ( Nl ) tR fN3, CELLS fCOMP, CHI, I DFfQtIER) 

DO  244  IC  =  1,3 

IF(Q.LT.CTST(IQ)  )NCHIN(N1,IQ)  =  NCHIN(N1,IQ)  +  1 

244  CONTINUE 

CALL    NKSHN0RM,R,N2,PDIF,IER) 

CC    245    IC=1,3 

IF(PDIFd)  .GT.STST(M,IQ)  INKSN(NltlQ)    =    MKSN(N1,I0)    ♦    1 

245  CONTINUE 

CALL    F0UTZ(NCRM,R,N2,Q) 

DO    246    10*1,3 

IF(0.GT.PNTST<N1,IC»)NFNN(N1,IC)    =    NFNN(N1,IQ)    +    1 

246  CONTINUE 
250  CONTINUE 
300    CONTINUE 

ISEED    =    CSEED 

WRITE(4»7)MEAN,VARI,KR,IS5P01,IS?EC 
WRITE(4,11)WT1,MEAN1,VAR1,WT2,MEAN2,VAR2 
11    FORMAT(//«    MIXED    NORMAL    (MIXNCPV)     =    » t F5  .2, ' *N ( • , F3. 1 , ■ , ■ , P3. 1, 
1    M    +«,F5.2,     «*N( «,F3.1,», •,F3.1,t )• //) 
WPITE(4,8)CHISCU,((NCHKN1.!Q)  . 

1  Nl=l,3)  ,10  =  1,31  ,<<NCHIN(N1,IC), 

2  N1=1,3),IQ=1,3) 
WRITE<4,8)K0LSNR,(<M<S<N1,!0),N1=1,3),IQ=1,3>, 

1    (  <NKSN<M,IQ>  ,N1  =  1,3)  ,10  =  1,3) 

WRITE(4,8)F0UTZL,((NFN<N1,I0),N1=1,3),I0=1,3) , ( ( NFNN ( Nl ,1 0) ,M1=1 
1    )  ,    10  =  1,3) 

7  FORMAT!//'    THE    VALLES    OF    MEAN     (U)     AND    VARIANCE     (V)    ARE1 ,2F12. 5// 
9    •    THE    RESULTS    OF    «,I5,«    REPL !CAT  ICNS,     STARTING    WITH    SEED*, 

1  I12//33X, 'ENDING    WITH    NEXT    SE=C',I12//) 

8  FORMAT(///20X,2A£//2GX, '20    PTS», 

2  6X,'30    PTS',6X,»50    PTS  V/10X  ,  » 10? ',  31 12/ •    CONTROL  *,  2X,  •    5*S 

3  3I12/10X,'    1*«,2I12///'    MIXNCPM'  ,2X,«10S» ,3112/'  VSN5X,'    5*' 

4  3112/'    N    <U,V)«,2X,«    1SS3I12) 

WRITE  (4, 9)  <(FNTST<M,IG),N1  =  1,3),IC=1,3) 

9  FORMAT(/«0  THIS    WAS    BASEC    ON    THE    FOLLLOWING    TEST    VALUES1, 
1    /(13X,3F12.5) ) 

GO    TO    ICO 
END 


SUBROUTINE    M IXNRM (ND,RAN,R ,DSEHD ,WT1 ,MEAN1,VAR1 ,MEAN2 , VAR2) 

IMPLICIT   REAL*8    (A-K,0-Z) 

DIMENSION    R(l)  ,RAN(1) 

REAL*4    RN(IOO) ,RNN,RNNN 

SVAR1     =    DSGRT(VAFl) 

SVAR2    =    DS0RT(VAR2) 

WT2    =    l.CO    -   WT1 
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XX   =   x*x 

TAN2    =     (PO    +    XX*<P1    «■    XX*P2M  /(PIBY4MQ0    +    XX*(Q1    +    XX*Q2)M 
RETURN 
200    TAN2    =    T4N(XARG)/XAFG 

RETURN 
END 

SUBROUTINE    FCUTZ     ( FCCF,XT ,NXT , FN ) 
C  THIS    SUBROUTINE    GENERATES    THE    STATISTIC    FOR   THE    FCUTZ    FN   TEST. 

C  INPUT    VARIABLES    ARE: 

C  PCDF    -    THE    CUMULATIVE    DISTRIBUTION    FUNCTION    AGAINST    WHICH    THE 

C  DEVIATES    ARE    BEING    TESTED.       CALLING    SEQUENCE    MUST    BE    OF 

C  THE    FCRM     •    CALL    PCDF<X,P)     •,    WHERE    X    IS    AN    INPUT    VALUE, 

C  AND    THE    VALUE    OF    THE    CUMULATIVE    DISTRIBUTION    FUNCTION 

C  IS    RETLRNSD    IN    P. 

C  P    MUST    BE    BETWEEN    0    AND    1. 

C  XT  THE    ARRAY    OF    DEVIATES,    IN    INCREASING    ORDER. 

C  NXT       -   THE    NUMBER    OF    DEVIATES    IN    THE    ARRAY    XT         (=    N    -    1) 

C 

C  THE    RETURNED    VALLE    IS    FN,    THE    VALUE    OF   THE    STATISTIC. 

C 

C  NXT    IS    PRESENTLY    LIMITED   TO    A    MAXIMUM    OF    50    BY    THE    DIMENSION    OF 

C  THE    VARIABLE    XTD. 

CIMENSION    XTC1) 

REAL*8    XTD(51> ,PN,FND 

N  =  NXT  4  1 

DO  200  1*1, NXT 

K  =  N  -  I 

CALL  PCDF(XT(K),P) 
200  XT0(K+1 I  =  P 

RN  -   l.DC/N 

XTD(l)  =  RN  -  XTC(2) 

DO  300  I =2, NXT 
300  XTD(I)  =  RN  -  XTCd+1)  +  XTD(I) 

XTD(N)  =  RN  -  l.CC  ♦  XTDfNl 

FND  =  0. 

CO  400  1*1, N 
400  FND  =  FNC  ♦  CMAX 1  (XTC < I ) , 0.00 ) 

FN  =  FND 

RETURN 

END 


39 


CALL    GGUBS(DSEED,2*ND,RN) 

CO    200    NN-1,ND 

RAN(NN)     -    RN(NN) 

CALL    MDNRIS<RN<NN>.RNN,IER> 

IF(RN(NN4ND).GT.WT2|G0  t0  150 

R(NN)  =  PNN*SVAR2  +  MEAN2 

GO  TO  200 
150  R(NN)  =  PNN*SVAR1  ♦  MEAN1 
200  CONTINUE 

RETURN 

END 


SUBROUTINE  NCRM(X,P) 

PEAL*8  MEAN.VARI ,SVARI 

COMMON/NOR PRM/M£AN,VARI,SVAP I 

T  =  (X  -  MEAN)/SVAPI 

P  =  .5*£PFC<-T*.7G71G68) 

RETURN 

ENO 


SUBROUTUE    UNIF<X,PI 

P    =    X 

RETURN 

END 


SUBROUTINE    FOUTZ    ( FCDF,XT,NXT , FN) 

C  THIS    SUBROUTINE    GENERATES    THE    STATISTIC    FOR    THE    FCUTZ    FN   TEST. 

C  INPUT    VARIABLES    APE: 

C  PCOF    -    THE    CUMULATIVE    DISTRIBUTION    FUNCTION    AGAINST    WHICH   THE 

C  DEVIATES    ARE    BEING    TESTED.       CALLIMG    SEQUENCE    MUST    RE 

C  THE    FCRV    •    CALL    PCC£<X,P)     •,    WHERE    X    IS    AM    INPUT    VALU 

C  AND    THE    VALUE    OF    THE    CUMULATIVE    DISTRIBUTICN    FUNCTION 

C  IS    RETURNED    IN    P. 

C  P    MUST    eE    EETWECN    0    AND    1. 

C  XT  THE    ARRAV    OF    DEVIATES,    IN    INCREASING    ORDER, 

C  NXT       -    THE    NUMBER    OF    DEVIATES    IN   THE    ARRAY    XT         (-    N   -    1) 

C 

C  THE    RETURNED   VALLE    IS    FN,    THE    VALUE    OF    THE    STATISTIC. 

C 

C  NXT    IS    PRESENTLY    UNITED   TO    A    MAXIMUM    OF    50    BY    THE    DIMENSION    OF 

C  THE    VARIABLE    XTD. 

DIMENSION    XT(1 ) 

REAL*8  XTD<51) ,RN,FND 

N  =  NXT  4  1 

DO  200  1*1, NXT 

K  =  N  -  I 

CALL    PCCF(XT(K),P) 
200    XTDCK+1)    =•    P 

RN    =    l.DO/N 

XTD(l)    =    RN   -    XTC(2) 

DO    300    1*2, NXT 
300    XTD(I)    =    RN    -    XT0U4-11    «•    XTD  C I  ) 

XTD(N)     =    RN    -    l.CO    ♦    XTC(N) 

FND    *    0. 

CO  400  1  =  1, N 
400  FND  =  FNC  ♦  0MAX1 (XTD < I ) ,Q  .DO  ) 

FN  =  FND 

RETURN 

END 
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c 

C  THIS  PROGRAM  GENERATES  PEARSON  TYPE  I  9R  II  RANDOM  DEVIATES. 

C  THE  PARAMETERS  ARE  CALCULATED  IN  TERMS  OF  Bl  ANO  B2  IN 

C  SUBROUTINE  PRM.   SUBROUTINE  ACINT1  IS  UScO  to  CALCULATE  THE 

C  CDF  TABLE.   GGUBS  IS  USED  TO  GENERATE  RANDOM  (CM  )  DEVIATES 

C  (UNIFORMLY  DISTRIBUTED,  AMD  THEN  SUBROUTINE  RANOP1  DOES 

C  AND  INVERSE  CALCLLATICN  TO  OBTAIN  THE  RANDOM 

C  PEARSON  CEVIATE.) 

C  A  SET  CF  50  RANCCN  CEVIATES  ARE  GENERATED,  THEN  THE  CHI 

C  SQUARED  TEST,  THE  KCLMOLGOROV-SMI PNOV  TEST,  AMC  THE  FOUTZ  <=N 

C  TEST  ARE  APPLIED,   THE  UNIFORM  DEVIATES  ARE  TESTED  AGAINST 

C  THE  HYPOTHESIS  THEY  CAME  FROM  A  UNIFORM  DISTRIBUTION  AS  A  CONTROL. 

C  THEN  THE  PEARSON  DEVIATES  ARE  Tg$TED  AGAINST  THE  HYPOTHESIS 

C  THAT  THEY  CAME  FPCM  A  NOPMAL  C TSTR TBUTION. 

C  THIS  TEST  IS  REPLICATED  A  NUMBER  OF  TIMES. 

C 

C  THE    TESTS    ARE    APPLIED   TO   THF    FIRST    20,    THE    NEXT    30,    AND   ALL    50. 

C 

C  THE   SET    OF    DEVIATES    IS    ALSO    TESTED   AGAINST    THE    NORMAL    DISTRIBUTION 

C  THIS    IS    REPEATED    FCR    THE   SAME    SETS. 

C 

C  THE  INPUT  VALUES  ARE  B1,B2,ISEC0,  AND  THE  NUMBER  OF  REPLICATIONS 


C 


IMPLICIT    REAL*8    (A-Z) 

INTEGER *4   T , IER ,N, I ,NP1,N1 ,N 2,N3 ,1  SEED , NR ,NSTOP,IDF , I Q , I  SEED! 

1  ,NCEL(3),NST(3> ,NSVP (3 ) ,NCH  I  (3,3 ) ,NCHIN ( 3, 3 ) , NKS< 3, 3 ) , 

2  NKSN( 3,3) ,NFNN(2,3)  ,NFN(3,3) ,NAFNN(3,3> ,NA-N ( 3 , 3 ) , NSP 
DIMENSION    X(2001),CDF(2001),A( 2001 ) , B( 200H ,RD ( 55) ,CHI SQU( 2) , 

1    K0LSMR(2)  ,P0UTZL(2)  ,ASFL(2)  , ^NTST (3 ,3 ) , ZALF( 3 ) , RN( 55 ) 

COMMON    /PPARM/C0,C1,C2,A1,A2,BC0,BC1,M1 ,M2 ,KINV,XL,XR ,MEAN,T 

COMMON    /CDF/X,CCF,A,E,NP1 

REAL*4   CELLS (50) ,CCNP (50)  ,  Q, CHI ,P ( 55 ) , PDIF(6 ) , OTST ( 3  ) , RNN ( 55  ) , 
1    STST(3,3) 

DATA    ST ST/. 26473,. 21 756,. 16959,. 29408,. 24170,. 18841,. 35 241, 
1    .28987, .22604/ 

CATA    NCEL/8, 10, 10/, NST/1, 21 , 1/ ,NSMP/20,30 , 50/ , QTST/. 1 ,. 05 , . 01/ 

1  ,CHISCU/«CHI    SQUA«,»RED    T=ST ■ /,KOLSMR/ 'KOLMOG-S ' , 'MIR    TEST1/ 

2  ,FOUTZL/«EMP    FOLT«,«Z    TEST* /  ,  ZA LCI , ZALC2/ .243069C0 , .36787900/ 

3  , ZALF / 1.2815500,1. 6 448 500,2.3 26  3 5C0/ 

OIMENSICN    FNEMP(3,3) 
DATA    FNEMP/. 42 714D0,. 41 903D0,. 4081600, 

2  .44965no,  .4355300, .4211600, 

3  .4865900,. 46579D0,.  4448700/ 

4  ,ASFL/»ASYM    FOU«,*TZ    T5STV 
EXTERNAL    NCRM,UNIF 

T   =    1 
100    PRINT    1 

READ   2,B1,B2 
IF(B1+B2.LE.0.D0)STCP 

2  F0RMAT(2E5.0) 
WRITE(4,14) 

14    FORMAT(lHl) 
1    F0RMAT(«4    Bl    AND    B2,     FORMAT ( 2E 5.0) ■ ) 

CALL    PRM(B1,B2,1.C0,ICR) 

IF(IER.NE.O)GO    TO    SCO 

A(l)    =    O.DO 

CDF(l)     =   0. 

CALL    INTSZ(N,X) 

NP1    -    N    ♦    1 

PRINT    3 

READ    5,  ISEED,NR 

ISEED1    =    ISESD 

DSEED    =    CFLOAT(ISEED) 
5     F0RMAT(I10,I5) 

3  FORMAT( ■♦INPUT    SEED    AND    NUMBER    OF    REPLICATIONS,    FORMAT( HO, 15 ) • ) 
CO    205    10=1,3 

DO    205    Nl*l,3 
NCHI(N1,I0)    '    0 
NCHIN(N1,I0)    -    0 
NKS(M1,IC)    =   0 
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NKSN(N1,I0>    =    0 

NFNtNl, IC)    =    0 

FNTST<N1,IG)    =    ZALCl*ZALF(IQ)/DS0RT<DFLOAT<NSMP(Nl)  +  in    ♦    ZALC2 

NAFNN<N1,IQ)    =    0 

NAFN(N1,IQ)    =    0 
205    NFNN<N1,IQ)    =    0 

DO    300    1=1, NR 

CALL    RANPDK50,RN,RC,0SEFD) 

DO    250    Nl=l,3 

N2    =    NST(Nl) 

N3    =   NSMP(Nl) 

CALL    VSPTAD<RD(N2),N3) 

CALL   VSRTAD(RN(N2),N3) 

DO    210    NSR=1,N3 

RNN(NSR)=    RNCNSR    ♦    N2-1) 
210    R(NSR)    =    RDCNSP+N2-1) 

IDF    =    0 

CALL    GFIT(UNIF,NCEL(N1),RNN,N2, CELLS, COMP, CHI ,IDF,Q,IER) 

CC    241    IC=1,3 

IFCO.LT.QTSTCIQ) )NCHI(N1,IQ)     =    NCHI(N1,IQ)    ♦    1 

241  CONTINUE 

CALL    NKSKUNIF,RNN,N3,P0IFtIEP.) 

DO    242    10    =    1,2 

IF(PDIFU)  .GT.STST<N1,IQ)  )NKS(N1,IQ>    =    NKS(N1,IQ)    +    1 

242  CONTINUE 

CALL    F0UTZ(UNIF,RNN,N3,Q) 

DO    243    IC=1,3 

I«=(0.GT.FNEMP(N1,IC)  )NFN(N1.IQ)    =    KFN(N1*I0)    +    1 

IFCQ.GT.FNTSTCNl,  101  »NAFN(N1,IQ1    =    NAFN(N1,IQ)    ♦    1 

243  CONTINUE 
IDF    =    0 

CALL    GFIT<N0RM,NCEL<Nl),R,N3,C=LLS,C0MPtCH!,IDe,Q,IER) 

DO    244    10    *    1,3 

IF(O.LT.QTSTUQ)  )NCHIN( Nl, ! Q )     =    NCHIN(  Nl  ,10)    ♦    1 

244  CONTINUE 

CALL    NKS1(N0PM,R,N3,PDIF,I=R) 

00    245    10=1,3 

IF(PDIF(1).GT.STSTCM,IQ))NKSN(M,IQ)    =    NKSN(N1,IQ)    +    1 

245  CONTINUE 

CALL    F0UTZ<N0RM,R,N3,Q) 

DO    246    10=1,3 

IF(0.GT.FNEMP(N1,IC) )NFNN(N1,TC)    =    NFNN(N1,I0I    +    1 

IF(G.GT.FNTST<N1,IC) )NAFNN(N1, 10)    =    NAFNN(N1,IQ)    +    1 

246  CONTINUE 
250  CONTINUE 
300    CONTINUE 

ISE5D    =    DSEED 

WRITE<4,7)B1,B2,A1,A2,M1,M2,NR,TSEED1,IS^0 
WRITE(4,8)CHISQU,((NCHI (N1,IC), 

1  N1=1,3),IQ=1,3),((NCHIN(N1,IQ), 

2  Nl=l,3) ,IQ=1,3) 

WPITE<4,£>KCLSNR,< (NKSCN1.IQI  ,  Nl=l ,3 ) , 10=1 ,3 ) , 
1    <<NKSN(N1,IQ),N1=1,3),IQ=1,3) 

WRITE(4,8)F0UTZL,((NFN(N1,IQ),N1  =  1,3),I0=1,3),  ( ( NFNN ( N1,I0 )  ,N1=] 
1    ),    10=1,3) 

WRITE(4,8)ASFL,(<NAFN(N1,IQ),N1=1,2),IQ=1,3), 
1    (<NAFNMN1,IC)  ,N1=1,3),I0=1,3) 

7  FORMAT<//»    THE    VALUES    OF    Bl    AND    B2    ARF» ,2F12 .5 /// 
A    •    VALUES    0*=    Al,    *2,    Ml,    AND    M2    ARE*  //4F12. 6/// 

9    •    THE    RESULTS    CF    »,I5,»    REPLICATICNS ,    STARTING   WITH    SEEDS 

1  I12//33X, 'ENDING    WITH   NFXT    SE=Df,I12//) 

8  F0RMAT<///2QX,2A8//2'3X,«20    PTSS 

2  6X,»30    PTS«,6X,«5G    PTS» //10X  , 'lO^1 ,3112/'    CONTROL* ,2X,  •    5?», 

3  3I12/10X,*    1**,3I12///*    P*AR$r)NN2X,  *10**,3I12/'         VS*,5X,*    5%  • 

4  3112/'     NORMAL*, 3X,*    1SS3I12I 
GO   TO    100 

900    PRINT    4,N,IFR,B1,E2,ERR0R 

4  FORMAT  (  *f\,IER,Bl,B2*  ,2110, 1P3C15.6  ) 
STOP 
END 
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SUBROUTINE  R*NPDl<NCfRANfRtOS5SD) 

IMPLICIT  R£AL*8  (A-^,C-Z) 
CIMENSICN  RIlliRANdl 
REAL*4  RN(10Q)tRNN 
CALL  GGUBS(DSEED,ND,RN) 
DO  200  NN=1,ND 
RAN(NN)  s  RN(NN) 
CALL  CDFINV(RAN(NN) ,RMN» 
R(NN)  =  PNN 
200  CONTINUE 
RETURN 
END 


SUBROUTINE  NORMX,F) 

P  =  .5*ERPC(-X*.7C71C68) 

RETURN 
END 


SUBROUTINE    UNIF(X,P) 

P    =    X 

RETURN 

END 

SUBROUTINE    PRM ( 81 t E2 tU2 , IER ) 

IMPLICIT   REAL*8    (A-Z) 

COMMON    /PPARM/C0,CltC2»Al,A2tBIGC0,BIGCl,Ml,M2»KINV?XL,XR,MEAM,T 

INTEGER*4    T,IEP 

IER    =    0 

IF(T.EG.2)B2    =    (6.    +    3.*Bl)/2. 

IF(T.E0.5)B1    =    4.*(2.*B2    -    6 . ) * (4.*B2-3 . )/ ( ( B2+3 . )**2+ 
1    12.*(4.*B2-2.  )) 

DEN    =    10.*B2    -    12.*ei    -    18. 

CO    =    <A.*B2   -    3.*ei)*U2/DEN 

CI    =    DSQPT(U2*B1)*<B2   +    3. I/DEN 

C2    =    (2.*B2    -    3.*E1    -   6.J/0EN 

Al    =    0, 

A2    =   0. 

BIGCO    =    0. 

BIGC1    =    C. 

n    -  0. 

M2    =   Q. 

GO   TO    (100,200,3GC,400,500,600,7CO), T 
100    CALL    RCCTS(C0,C1»C2,A1,42,  IER) 

IF(IER.NE.O)RETUPN 

Ml    =    (CI    ♦    A1)/C2/(A2   -   Al) 

M2    =    -(CI    +    A2)/C2/(A2    -   Al ) 

KINV    =     (A2   -    A1I**(M+    M2    +    1.  )*DGAMMA(MH-1 .  )  *DGANMA  (  M2+1  .  ) 
1    /DGAMMA(M1    ♦    M2    4    2.) 

XL    =    Al 

XP    =   A  2 

GO  TO  750 
200  GO  TO  100 
300    Ml    =    (C0/C1   -    Cll/Cl 

C2    =   0. 

KINV    =    DGANMA( N1+1.I*DEXP(C0/C1**2)*DABS(C1)**(2.*M1+1.) 

XL    =    -C0/C1 

XR    =    1.C50 

GO  TO  7  5C 
400    BIGCO    =    CO   -    Cl**2/C2/4. 

BIGC1    =    C1/C2/2. 

Ml    =    (CI   -    3IGC1)/0SQRT(2.D0)/9IGC0 

M2    =    DSQRT(BIGC0/C2) 

KINV    =    l.DO 

XL    -    -1.D50 

XR    =    1.D50 

GO  TO  7  5C 
500    BIGC1    =    C1/2./C2 

43 


Ml    =    (CI    -    BIGCU/C2 

IF(C2.GT.1.DC)G0    TO    800 

KINV    =    CABS(M1I**(1./C2   -    1. ) /OGAMMA ( 1. /C2-1.  ) 

XL    =   -BIGC1 

XP    =    1.D50 

GO   TO   750 
600   CALL    R0CTS(C0,C1,C2,A1,A2,IER) 

Ml    =   -(CI    ♦   A1)/C2/(A2   -   Al> 

M2    '    (CI    ♦    A2)/C2/(A2    -   Al  ) 

IF(M2.GE.-1.D0.0R.P1    +    M2 . G5.0  .00 IGQ   TO    800 

KINV    =     (A2   -    A1)**(M1    +    M2    +    1. )*DGAMMA( M2+1. ) *DGAMMA (-M1-M2-1 . ) 
1    /DGAMM(-Ml) 

XL    =   A  2 

XR    =    1.050 

GO    TO    750 
700    KINV=C0**(-.5D0/C2)*DSQRT(C0/C2)*DGAMMA(.5D0)*DGAMMA(  .5D0/C2-.5C 
1    /0GAMMAU500/C2J 

CI    =  0. 

XL    =    -1.050 

XP    =    1.C50 
750    MEAN    =     (4.D0*B2    -    3 .D0*B1 ) /3.D0/ ( B2    -    Bl    -   1 . DO l*U2*PDF(0 .00 ) 

RETURN 
800    IER    =    1 

RETURN 

END 

SUBROUTINE    ROOTS (CO, CI, C2, A  1, A?, IER) 

IMPLICIT    REAL*8     (A-Z) 

INTEGER**    IER 

IER    =    0 

IF(C2.EC.0.D0IGC    TC    500 

DIS    =    Cl**2   -    4.CC*C0*C2 

IF(DIS.LT.O.DO)GC    TC    600 

SDIS    =    CSQRT(QIS) 

DNUM    =    -CI    -   SDIS 

IF(C1.LT.0.D0)CNUN    =    -CI    ♦    SCIS 

X2    =   DNLM/2.DO/C2 

XI    =   C0/C2/X2 

Al    =  DMINHX1.X2I 

A2    =    DMAX1(X1,X2) 

RETURN 
500    A2    =    1.D75 

Al    =    -C0/C1 

RETURN 
600    IER    =    -1 

RETURN 

END 

FUNCTION    PDF(X) 
C 

C  THIS    FUNCTION    EVALUATESTHE    PEARSON    DISTRIBUTION 

C  FOR    A    GIVEN    X,    THE    PARAMETERS    HAVING    PREVIOUSLY    BEEN 

C  CALCULATED    IN    SLBRCLTINE    PRM. 

C 

IMPLICIT    REAL*8    (A-Z) 

INTEGER**    ITY 

COMMON/ PPARM/    CO ,C1 ,C2, Alt  A 2 , B !GC0,BIGC1, M1,M2 ,KINV,  XI ,X2 ,MEAN  ,1 

GO    TO    (100,200,  300, 400 ,500 ,600 ,790 )♦ !TY 
100    IF(X.LE.A1.0R.X.GE.A2)G0   TO    140 

PDF    =    (X    -   A1)**N1*(A2    -    X)**M2/KINV 

RETURN 
140    IF(X.LE. ALAND. Ml. LT.O.DOJGO    TO    15C 

IF(X.GE.A2.AND.M2.LT.0.C0)G0    TO    15C 

PDF    s    l.C-25 

RETURN 
150    PDF    =    1.C25 

RETURN 
200    GO   TO    IOC 
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300    PDF    =    (CO    +    C1*X)**M*DEXP(-X/C1)/KINV 

RFTURN 
400    PDF   =    (BIGCO    ♦    C2*  (X+BIGC  1 )  **2  )**(  -1.  /C2)  * 
1    DEXP(-M*DATAM  (X    ♦    PIGC1 )  /  f*2  )  )  /K  INV 

RFTURN 
500    PDF    =    (X    +    BIGC1I**(-1./C2)*D*=XP(M1/(X    ♦    BIGC1U/KINV 

RETURN 
600    PDF    =     (X    -    A1)**M1*(X   -   A2)**M2/KINV 

RETURN 
700    PDF    =    (CC    ♦   C2*X**2)**(-.5D0/C2)/K!NV 

RETURN 

END 


SUBROUTUE    INT5Z(N,X) 

IMPLICIT    REAL*8    (A-Z) 

COMMON    /FPARM/C0,Cl,C2,Al:A2,eCQ,BCl,Ml,M2,KINV,XL,XR,MEAN,T 

COMMON /CCF/XDUM( 2001) ,CDF ( 2001 ) t A ( 2001 ) ,6(2001 ),NP1 

CIMENSICN    X(l) 

INTSGER*4    T,N,NN,INT,IEND,NP1,IER 

X<1)    =    Al 

B(l)    =    PCF(Al) 

CDF(l)     -   O.DO 

INT    =    1 

NN    ■  1 

DXM    =    (A2    -   All/1. C7 

CX    =    DXM*1.D3 

GO    TO    160 
150    DX    =    DX*(1.D0   +DSGRT(.475D-4/=REST))/2. 

INT    =    1 
160    X(NN+1)     =    X(NM    4    CX 

IEND    =    1 

IF(X(NN4l) .I.T.A2  )G0    TO    170 

X(NN+>1)     -    A2 

DX    =    X(NN+1)    -    X(NN) 

IEND    =    2 
170    B(NN+1)     =    PDF(X(KK+1II 

ESTINT    -    (B(NN-H)    ♦    B  (NNM  *D  X/2.D0 

IF(ESTINT.LT.1,C-8)G0   TO   301 

300  DCDF    =    ADINTK  X(NM  ,X(NN+i)  ,EQRCR,I5R) 
IF(IER.GT.100)G0    TC    900 
IF(DCDF.GT..025D0)GC    TO   309 

EREST    =    (l.D0/B(NN+l)    -    1. DO/B (NN)  )*DCDF/8  .DO 
A(NN)    -    EREST 
EPEST    =    CABS(ERSST) 
GO   TO    (202,308),  INT 

301  DCDF    =    ESTINT 

EPFST    *    5.27777777170-6 
A(NN)    =    EREST 
GC    TO    308 

302  IF(EREST.LE..5C-4)GG    TH    308 

CX    -    DX*(DSQRT(.25C-4/EREST)    ♦    l.DC)*.5D0 
IF(DX.GT.DXM)GC    TC    304 
DX    =    DXM 

X(NN+1)     =    X(NN)    ♦    CX 
INT    =    2 
GO   TO    17C 
304    X(NM-l)     =    X(NN)    *CX 
IEND    =    1 
GO   TO    170 

308  CDF(NN-H)    =    CDF(NN)    +    DCDF 
NN    =    NN    +    1 
IF(NN.GT.2000)G0    TO    910 

GO    TO    (150,310)  ,IEND 

309  DX    =    DX*.0125DG/CCDF 
GO   TO    160 

310  N    =   NN    -    1 
NP1    =    NN 

DO   400    Nf\  =  l,N 
400    CDF(NN+1)    =    CDF(NN+1)/CDF(N+1) 
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2    FORMAT <•    CDF(NN)    =     »,F20.15) 
RETURN 
900    PRINT    l,NNtX(NNI ,)(NN+1I ,EFRCR 
STOP 
1    FORMATC    I  ,XI  ,XI-H,ERPCR»  ,  I  5  , 1 P3C15.6 ) 
910    PRINT    4,X(NN) 

4    FCRMAT(J       RUN    CUT    OF    SPACE    AT    X=«,F12.5) 
STOP 
ENO 


SUBROUTINE    CDFINV(CP) 

R£AL*8    X,C0F,A,B,CC,CCC,ADINT1,XC,E 

COMMON/CDF /X( 2001) tCCF (2001 )  ,  A (2001 > , B(2001 ) , NP1 

CALL    TABLOC(CDF,CeLE(C),NPl,I) 

XC    =    X(I)    +    (C-CCF(I) )/(CDF(I*l)    -    CDF(I) )*(X( 1+1)    -    X(!)) 

IF(DABS(A(in.LE..5D-4)G0    TO    410 
300    CO   400    J=l,10 

DC    =   ADINTHXdl  ,XC,E,IEP) 

DOC    =    (CDF(I)    ♦    DC    -    C)/POF(XC) 

XC    =   XC    -    DOC 

IF(DABS(CDC).LT..5C-4)G0    TO    410 

IF(XC.LE.X(T))GO    TO    350 

IF(XC.GE.X<t*l  )  )GC    TO   360 

GO    TO    4CC 
350    XC   =    X(I)    ♦    (C    -    COF(I) )/DC*(XC    -    XII I  I 

360    XC    =    X??+ll    +    (C    -    CDF(I+1) )/(COF(II    +    DC    -    CDF( 1+1 ) )* (X ( 1+1 )-X 

400    CONTINUE 

URITE(6,1)C,XC 
1    FORMATS    10    ITERATICNS    IN    CDFINV,    C    AND    XC    =»,2E12.4) 
410    P    *    XC 

RETURN 

END 


SUBROUTINE    TABLOC ( XT, X, N ,NL ) 

REAL*8    XT(1),X 

NT    =   ALGG(FL0AT(M)/.301    ♦    1. 

NU    =    M 

NL    =    1 

DO    200    1=1, NT 

NG  *    (NU    ♦   ND/2 

IF(X.GE.XT(NG) )GC    TC    103 

NU    '    NG 

GO   TO    200 
100    NL    =    NG 
200   CONTINUE 

RETURN 

END 

FUNCTION   ADINT1(X1.X2,ERR0P,IEP> 
C  THIS    SUBROUTINE    USES    DCADR5    TO   OBTAIN    THE    INTEGRAL    OF    A 

C  PEARS3N    TYPE    I    OR    II    DISTRIBUTION    FUNCTION,       SUBTRACTING 

C  OUT   THE    SINGULARITY    IS    USED. 

C 

IMPLICIT  REAL*8  IA-Z) 

INTEGER**    T,IER 

EXTERNAL    F1,F2,PCF 

COMMON/PP4RM/CO,C1,C2,A1,A2,BCO,BC1,M1,M2,KINV,XL,XR,MEAN,T 

IF(Xl.GT.(XL+XR)/2.CQjG0   TO    2C0 

IF(M1.GT.0.D0)GC    TC    300 

PR    =   DCADRE(F1, XI, X2  ,1 . C-6 ,0.C0, ERROR, IER) 

ADINT1    =    PR    ♦    <A2-A11**M2*UX2-A1I**(M1«-1.D0)-<X1-A1)**(M14-1.D0 
1    /(Ml    +    l.DO)/KINV 

RETURN 
200    IF(M2.GT.0.D0)G0    TC    300 

PR    =   DCACRE(F2, XI,  X2, 1.0-6,0. CO, ERFOR,IERI 

A0INT1    =    PR-(A2    -    A1)**M1*UA2-X2)**(M2+1.00)-<A2-X1I**(M2«-1.DO 
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1    /<M2+1.C0)/KINV 

RETURN 
300    ADINT1    =    0CADRE(PDF,X1,X2, 1 .D-6,0. CO, ERROR ,IEP ) 
RETURN 
END 

FUNCTION    FKX) 
IMPLICIT    RSAL*fi    (A-Z  ) 

COMMON    /FPARM/C0,Cl,C2,Al,A2,eC9,BCl,Ml,M2,KINV,XL,XR,MEAN,T 
IFCX.LE.A1IG0    TO    20C 

Fl    =•    TX-An**Ml*UA2    -    X)**M2    -    (A2   -   Al) **M2) /KINV 
RETURN 
200    Fl    =    O.DC 
RETURN 
ENO 

FUNCTION    F2(X) 
IMPLICIT    REAL*8(A-Z) 

CGMMON/ F F ARM/CO, CI,  C2,  A1,A2,8C0,8C1, Ml, M2, KINV, XL, XR, MEAN, T 
IF(X.GE.A2)G0    TC    2C0 

F2    =    <A2-    X)**M2*(<X-A1)**M1    -    (A2-A1) **M1) /KINV 
RETURN 
200    F2    =    O.DC 
RETURN 
END 

SUBROUTINE    FCUTZ    ( PCDP ,XT,NXT  ,  FN  ) 
C  THIS    SUBROUTINE    GENERATES    THE    STATISTIC    PHR    ^HS    PCUTZ    FN   TEST. 

C  INPUT    VARIABLES    ARE: 

C  PCDF    -    THE    CUMULATIVE    DISTRIBUTION    FUNCTION    AGAINST    WHICH   THS 

C  DEVIATES    ARE    BEING    tcSTcD#       CALLING    SEQUENCE    MUST    Bc    OP 

C  THE    FCPM    ■    CALL    PCCP(X,P»    •  ,    WHERE    X    IS   AN    INPUT    VALUE, 

C  AND    THE    VALUE    OF    THE    CUMULATIVE    DISTRIBUTION    FUNCTION 

C  IS    RETUPNEC    IN    P. 

C  P    MLST    EE    BETWEEN    0    AND    1. 

C  XT  THE    ARRAY    OF    DEVIATES,    IN    INCREASING    ORDER. 

C  NXT       -   THE    NUPEER    CP    DEVIATES    IN    THE    ARRAY    XT         (=    N   -    1) 

C 

C  THE    RETURNED    VALUE    IS    PN,    THE    VALUE    OF    THE    STATISTIC. 

C 

C  NXT    IS    PRESENTLY    LIMTED    TO    A    MAXIMUM    OF    50    BY    THE    DIMENSION    OF 

C  THE    VARIABLE    XTO. 

DIMENSICN    XT(1 ) 

REAL*8  XTDC51) ,RN,FND 

N  =  NXT  ♦  1 

DO  200  1=1  ,NXT 

K  -   N  -  I 

CALL    PCDF(XT(K),P) 
200    XTDCK+1)     -    P 

RN    =    l.DC/N 

XTD(l)     »    RN    -    XTCC2) 

DO    300    I =2, NXT 
300    XTD(I)    *    RN   -    XTC(I-H)    +    XTD  ( I! 

XTD<N»     »    RN    -    l.CO    ♦    XTO(N) 

FMD    =    0. 

CO  400  1=1, N 
400  FND  =  FNC  ♦  CMAXKXTC  ( I )  ,3  .DO  ) 

FN  =  FND 

RETURN 

ENO 
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APPENDIX   2 

CHI     SQLUREC    TEST 

20    PTS  30    PTS  50    FTS 


CCN7R0L 

1C? 
52 
1? 

410 

2  40 

48 

429 

240 
55 

461 

229 

56 

RSSF 
VS 
N    (0,1) 

10? 
5% 
12 

1784 

1362 

522 

2908 
2384 
1410 

4125 
3645 
2673 

KOLMOG-SMIR    TEST 

2C    PTS  30    PTS  5C    PTS 

102                         491  495  463 

CCNTFQL         52                        253  253  224 

12                           49  58  46 

RSSF            1C?                      1636  2202  3316 

VS                 52                       1008  1401  2410 

N    (0,1)         1?                         270  466  9C5 


FQUTZ    FN    TEST 


2C    PTS 

30    PTS 

50    PTS 

CONTROL 

10? 
5% 
1? 

525 
27C 

4C 

481 

250 

35 

488 

253 

62 

RSSF 
VS 
N    (0,1) 

102 
52 
12 

36  03 
2999 
1894 

4010 
3515 
2463 

4597 
4279 
3468 

•CCNTRCl'     IS    TFE    TEST    OF    LNIF0RM(0,1)    VS    UNIFCRM(C.l) 

•RSSF'     IS    THE    F4ND0M    STAEILIZEO    STANDARD    FORM 
XLPHA    =    1.00,    eETA     =    CO 


NUMBER    OF    HYPC7HESIS    REJECTIONS     IN    5000    REPLICATIONS 
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CHI     SQLAR  ED    TEST 
20    PTS  3G    ojS 


5C    PTS 


CCNTFCL 

10? 

5$ 
1? 

39C 

220 

49 

421 

219 

57 

498 

242 

54 

RSSP 
VS 
N    (  C,l  ) 

1C? 
5? 
1? 

1460 

1084 

401 

2407 
1889 
1011 

3616 
3039 
1967 

CCNTPCl 


CCNTFGL 


10? 
5* 
1? 


RSSF  10? 

VS  5? 

N    <0,1)  1? 


10? 
5? 
1* 


RSSF  10? 

VS  5? 

N    (0,1)  1? 


KCLMOG-SMIR     TEST 
20    PTS  3C    PTS 


488 

260 

64 

46C 

244 

50 

1412 
792 
222 

1782 

1083 
324 

FOUTZ    FN 

TFST 

20    PTS 

3C    PTS 

503 

2  56 

60 

494 

257 

43 

28  6  3 
2145 
11  C2 

3272 
2678 
1575 

50  PTS 

500 

231 

38 


2654 

1760 

57« 


50  PTS 

5C1 

241 

42 


3963 
3422 
2244 


•CCNTRCL1      IS    THE    TEST    CF    UNIFQRM(0,1)     VS    UNIFORN(Cl) 

•RSSF»     IS    THE    FANDOM    STAEILIZEC    STANDARD    FORM 
4LFH/S    =    1.30,     EET4    =    0,0 


NIMEER    OF    HYPOTHESIS    REJECTIONS     IN    5000    REPLICATIONS 
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CHI    SQU/RFD    TFST 

20    PTS  3C    PTS  5C    FT< 

10?                         3  95  29  5  5C5 

CCfvTFCL         5%                        221  199  231 

13                            39  40  45 

RSSF            1C?                       1441  2116  3266 

V<                 53                       1054  1598  2692 

N    (C,l)          1?                         390  817  1663 


KOLMOG-SMR     tcST 

20    PTS  30    PTS  50    PTS 


CCNTRCL 

10? 
5% 
1% 

493 

245 

5C 

472 

215 

47 

5C2 
265 

52 

PSSF 
VS 
N    (0,1) 

103 
53 
13 

1321 

780 
210 

FQUTZ    Fiv 

1635 
991 
294 

TEST 

2362 

15C8 

521 

20    PTS 

3C    PTS 

50    FT 

CCNTRCL 

1C? 

53 
13 

468 

188 

27 

449 

222 

51 

496 

237 

40 

PSSF 

VS 

N    (0*1) 

10? 
5? 
1% 

2479 

18C4 
852 

2678 
2059 
1080 

33C6 
27C5 
1572 

•CONTROL1     IS    THE    TEST    CF    UNIFQRM(0,1)    VS    UNIF0RM<C,1) 

•RSSF*     IS    THE    RANDOM   STABILIZED    STANDARD    FORM 
ALPH4    ~    1.60,     EETA    -    0.0 


MMEEP    OF    HYPCTHESIS    REJECTIONS     IN    5000    REPLICATICNS 


50 


CHI    SQUARED    TE5T 

2C    °TS  30    PTS  50    PTS 

10?          389  4C5  491 

CONTROL    5?         230  221  235 

1?           43  47  4£ 

PSSF            10?                       1326  1933  3046 

VS                 5%                         932  1469  2428 

N     (0,1)          1%                         32C  688  1440 


KGL^OG-Si^IR    TEST 

20    PTS  3C    PTS  50    PTS 


CCN7P0L 

1C? 

53 
IS 

5  27 

275 

59 

474 
236 

38 

489 
262 

54 

RSSF 
VS 
b    (0,1) 

10? 
13 

1226 

760 
242 

1521 
9  29 
271 

2184 

1389 

474 

FOUTZ  FN  TEST 

20  PTS  30  PTS  50  PTS 

1CS         531  505  516 

CCNTPGL    5%                         265  263  262 

1%                            53  49  54 

RSSF            1C?                       2046  2274  2716 

VS                 5%                      1396  1646  2061 

N    (C,l)         13                        581  78£  1092 


'CCNTPCL*    IS    THE    TEST    OF    UNIF0PM(0,1)     VS    UNIF0PM(0,1) 

•RSSF*     IS    THE    RANDOM    STAEILIZEC    STANDARD    FOR f 
ALPHA    *    1.90,     BETA    =    0.0 


MJNBEF    OF    HYPCTHESIS    REJECTIONS     IN    5000   REPLICATIONS 


51 


CHI     SQUARED    TEST 

20    PTS  3C    PTS  50    PTS 

iCS                        400  457  4*56 

CONTROL         5%                         234-  242  246 

1%                           40  51  57 

PSSF            103                       1728  2935  4056 

VS                 5%                       1316  2435  3629 

N    (0,1)         1%                        555  1454  2672 


KOLMQG-SMR    TEST 

20    PTS  3C    PTS  50    PTS 


CONTROL 

1C3 

5% 
1% 

504 

252 

68 

474 

235 

52 

477 

234 
47 

FSSF 
VS 
N    (0,1) 

10? 

5% 
1% 

1565 
966 
31C 

FOUTZ    FN 

2222 

1410 

481 

TEST 

3221 

2372 

901 

20    PTS 

3C    PTS 

50    FT 

CONTROL 

1? 

493 

238 

49 

522 

255 
53 

467 

250 

46 

RSSF 
VS 
N    (0,1) 

103 

5% 
1% 

3502 
2859 
1753 

4014 
2522 
25^7 

4596 
42S9 
3487 

•CCNTFOL'     IS    THE    TEST    0*=    UNIF0RM(0,1)    VS    UNIF0RM(0,1) 

'P-SSF*     IS    THE    RANDOM    STAEILIZEC   STANDARD    FCRM 
ALPHA    =    I. CO,     eETA    =    0. 25 


K'NBEF    GF    HYPCTHESIS    REJECTIONS    IN    5000    REPLICATIONS 


52 


CHI     SCUARED    TEST 

2C    PTS  30    PTS  50    PTS 

10X          372  416  465 

CONTROL    5%                        225  232  204 

IS           42  56  45 

RSSF     103         1521  2451  3593 

\S       5%                      1128  1933  31C6 

N  (0 ,1 )    1*          454  1057  2046 


KOLMOG-SMIR     TEST 

20    FTS  3C    PTS  50    PTS 

10?                        486  512  483 

CCNTRCL         51                       242  272  25E 

135                           48  62  45 

RSSF            1C*                      1430  1849  27C1 

VS                 5%                        890  1171  1882 

N    (Otl)          1*                        268  371  67< 


fCUTZ    FN    TEST 

20    PTS  30    PTS  50    PTS 

1035                        461  465  459 

CONTROL         5%                       232  252  228 

1%                           37  53  47 

PSSF            1055                      2868  3284  4026 

VS                  53                      2236  2684  3486 

N     (0,1)          1%                      1174  1545  2329 


•CCNTRCL'     IS    THE    TfST    OF    UNI  FORM  (  0  ♦  1)    VS    UNIFORNCGrl) 

•RSSF*     IS    TVE    RANDOM    STABILIZED    STANDARD    FORM 
ALPH/S    =    1.  30t     EETA    =    C.25 


NLMeEP    OF    FYPC1HESTS    REJECTIONS     IN    5000    REPLICATIONS 


53 


CHI    SQUARED    TEST 
20    PTS  3C    PTS 


50    PTS 


CONTROL 

1C3 
1% 

395 

2  23 

31 

390 

212 

54 

467 

2  42 

53 

PSSF 
VS 
N    (CD 

10% 
5% 
12 

14  16 

1C39 

357 

2100 

16  10 

786 

3277 

27C5 
16C7 

CCNTFOL 


10? 
5% 
1% 


KOLMGG-SM1R    TEST 
20    PTS  30    PTS 


514 

260 

50 


518 

282 

61 


5C    PTS 

504 

278 

65 


RSSF  1C3 

VS  5% 

N    (  CI  I         1? 


1345 
807 
239 


1661 

1C17 

331 


24C2 

1616 

539 


FOUTZ    FN    TEST 


2C    PTS 

30    PTS 

50    PTS 

CCNTPCL 

10  % 

5% 
1% 

49  2 

248 

45 

5C2 

260 

35 

495 
24C 

C  1 

PSSF 
VS 
N    (0,1) 

10? 
5? 
1% 

23  83 

1698 

822 

2650 
2023 
1027 

3260 
26  11 
1413 

'CCNTRCL*     IS    THE    TEST    CF    UN!FCRM(0,1)    VS    UNIFORM(0,1) 

•RSSF'     IS    THE    RANDOM    STAEILIZEC    STANDARD    FORM 
^LPH/    =    1.60,    eETA     -    C.25 


NUMeEP    OF    HYPCTHESIS    REJECTIONS     IN    5000    REPLICATICNS 


54 


CHI     SQUARED    TEST 

20    PTS  3C    PTS  50    PTS 

1C^!          402  441  497 

CCNTROL    5*         222  230  237 

IS                           43  45  46 

RSSF            10?                      1289  1916  3  044 

VS                 5%                        933  1458  2409 

N    (0,1)         1%                        319  652  1326 


KCLVCG-SMI  R    test 

20    PTS  3C    PTS  50    PTS 

103                          488  469  487 

5%                        234  235  242 

i%                           40  52  52 

PSSF            10?                      1239  1471  2111 

VS                 5%                        713  904  1357 

N     (0,1)          1?                        223  269  445 


488 
234 

40 

46  9 

235 

52 

1239 

713 
223 

1471 

904 
269 

FOUTZ    FN 

TEST 

20    PTS 

3C    PTS 

484 

215 

56 

470 

246 

50 

20  0  3 

1358 

552 

2289 

1688 

770 

50    PTS 

1C  -  470  494 

CONTROL         5?  215  246  242 

13  56  50  49 

RSSF           103?  ,                 2003  2239  2773 

VS                 5%  1358  1688  2098 

h    (0,1)         1*  552  770  10C7 


•CCM7RCL'     IS    THE    TEST    OF    UMFCRM(O-l)    VS    UNI  FORM  (0, 1  ) 

•RSSF»     15    THE    F-ANDOM    STAEIUZEC    STANDARD    FORM 
ALPHA    =    1,90,     EETA    =    0.25 


M'MEEP    OF    HYPCTFESIS    REJECTIONS    IN    5000    R  EP  LICAT  IC  f^S 


55 


CCMPCl 


10? 
5% 
1* 


CHI    SOU4R  ED    TEST 
20    PTS  30    PTS 


3  83 
213 

45 


427 

248 

56 


50     FTS 

49<5 

253 

37 


RSSF 
VS 
IS    (0,1) 


1035 

5% 
1% 


1732 

13C3 

539 


2905 
2395 
1432 


4113 
3668 
2660 


KCtVCG-SMIR    TEST 
20    PTS  30    PTS 


50    PTS 


CCNTFCL 

10* 

5? 
1% 

5C7 
2  45 

38 

49  2 

233 

47 

4£4 

232 

37 

RSSF 
N  (GU) 

1C3 

5% 
1? 

1663 
1011 

289 

FOUTZ  FN 

2188 

1430 

455 

TEST 

3292 

2414 

93C 

2C  PTS 

30  PTS 

50  PT 

CCNTPCl 

102 

5% 
1% 

507 

2  59 

53 

492 

237 

58 

534 
271 

50 

RSSF 
VS 
N  (0,1) 

10* 
5* 
1* 

34  92 
2399 
1822 

4011 
3512 
2489 

4571 
4297 

346C 

'CCNTRCL1     IS    THE    TEST    D=    CNI FORM (0  , 1 )    VS    UNIFCRM(C,1) 

•RSSF'     IS    THE     F4ND0M    STAEILIZED    STANDARO    FORM 
^LPH*    =    l.CO,    BETA     =    0,50 


NUNeER    OF    HYPOTHECS    REJECTIONS     IN    5000    REPLICATICNS 


56 


CHI     SQLAREO    TEST 

20    PrS  3C    PTS  50    PTS 


CONTROL 

1C? 
5% 
1% 

382 

213 

45 

406 

19  5 

38 

442 
2C7 

3€ 

RSSF 
VS 
N    (0,1) 

10? 
5? 
1? 

1665 

1266 

509 

2552 
2  03  8 
1178 

36S5 
3222 
22C2 

KCLMCG-SMIR    TEST 

20    PTS  30    PTS  50    PTS 

10*                         501  503  486 

CONTROL         5$                        248  256  253 

1%                            45  55  57 

PSSF             10$                       16CC  2092  3025 

VS                 5%                       1026  1425  2208 

tv     (0,1)          It                         342  532  951 


501 

248 
45 

5  03 

256 

55 

16CC 

1026 

342 

2092 

1425 

532 

FOUTZ    FN 

TEST 

20    PTS 

3C    PTS 

5  06 

246 

44 

467 

236 
46 

28  16 
2168 
1180 

2278 
2680 
1546 

50    PTS 

1C  5  06  471 

CrNTROL         5%  246  236  228 

1?  44  46  43 

RSSF           10?  2816  2278  392C 

VS                  5?  2168  2680  3382 

N     (0,1)          1?  1180  1546  227C 


•CCNTRCL*     IS    THE    TEST    OF    UMFCRM(0,1)    VS    UNIF0RM(0,1) 

•RSSF»     IS    IFF    RANDOM    STAeiLIZEC    STANDARD    FORM 
ALPHA    =    1.30,     EETA    =    0. 50 


M.MBE  F    CF    HYPCTFESIS    REJECTIONS    IN    5000    REPLICATIONS 


57 


CHI    SQUARED    TEST 


20    FTS 

3C    PTS 

50    PTS 

CONTROL 

1C3 
5? 

1% 

413 
251 

36 

440 
222 

48 

470 

2C5 

41 

RSSF 
VS 

N    (0,1) 

1C3 
5* 
1? 

14  51 

1063 

366 

2269 

1794 

948 

3435 
2841 
1626 

KCLNCG-SMI R    TEST 

20    PTS  30    PTS  50    PTS 

10?          478  518  513 

CONTROL    5%                         252  278  277 

1%                            49  65  6C 

PSSF            10?                       1435  1853  2636 

VS                 5%                         888  1229  1881 

N     (0,1)          15g                        267  444  738 


478 

252 

49 

518 

278 

65 

1435 
888 
2  67 

1853 

1229 

444 

FOUTZ    FN 

TEST 

20    PTS 

3C    PTS 

493 

241 

47 

463 

22C 

47 

2314 

16  86 

761 

2643 

202  8 
1C52 

50    PTS 

1C?  493  463  488 

CCNTFOL         5%  241  22C  227 

1%  47  47  40 

RSSf            101  2314  2643  3268 

VS                  5%  1686  2028  2594 

N     (0,1)          1?  761  1C52  1449 


•CCN7PCL*     IS    T>E    TEST    OF    UMFCRM(0,1)    VS    J  NI  FORM  (0,  1  ) 

•RSSF*     IS    THE    RANDOM    STAETLIZEC    STANDARD    FORM 
ALPHA    =    1.60,     EETA    =    0*50 


NLMEER    OF    HYPOTf-ESIS    REJECTICNS    IN    5000    REPLICATICNS 


58 


CHI    5QU4PED    TEST 

20    PTS  20    PTS  5C    PTS 


CCNTFQL 

10* 

52 
1% 

419 

218 

38 

411 

230 

54 

466 

224 

52 

RSSF 
\S 
N    (CD 

1C3 

5% 
13 

1299 
916 
297 

1934 

1438 

680 

3C42 
24C7 
1348 

KOLMOG-SWR     TEST 

20    FTS  30    PTS  50    PTS 

10%          5  C2  494  48  2 

CONTPCL    5?          26C  255  24g 

IS           47  58  62 

PSSf            1C?                       1239  1559  214E 

VS                 5%                        738  961  1357 

N    (C.l)          1?                        214  296  475 


5C2 

26C 

47 

494 

255 

58 

1239 
738 
214 

1559 
961 
296 

FCUTZ    FN 

TEST 

2  0    PTS 

30    PTS 

473 
2  59 

479 
229 
53 

2018 

1384 

591 

2310 

1629 

749 

50    PTS 

103  508 

CONTROL         51  259  229  240 


PSSF  10?  2018  2310  2708 

VS  5%  1384  1629  2043 

N     (0,1)  1?  591  749  1023 


•CCNTRCL'     IS    THE    TEST    OF    UNIF0RM(0,1>    VS    UNIFORMED 

•RSSF*     IS    TFE     FANDOM    STAEILIZEQ    STANDARD    FORM 
ALPH*    =    i,90,    eETA    =    0.50 


MMBEF    OF    FVPCHESIS    REJECTIONS     IN    5000    REPLICATIONS 


59 


CHI    SQUARED    TEST 

20    PTS  3C    PTS  50    PTS 

1C?          410  445  436 

CONTROL    1%                         229  236  221 

1%                          30  49  46 

RSSF            10?                      1831  2913  41C8 

VS                  5%                       13*8  T391  365S 

*     (C,li         1%                         566  1438  2690 


KOLMOG-SMIR    T  EST 

2C    PTS  30    PTS  5C    PTS 


CCNTFQL 

10S 
5? 
138 

524 

250 

46 

505 

259 

53 

527 

25  1 

5<5 

PSSf 
VS 
N    (0,1) 

icsg 

53 
1? 

1627 
10  51 

3  04 

FCUTZ    FN 

2196 

142  1 

500 

TEST 

33C4 

2387 

931 

2C    PTS 

30    PTS 

50    PT 

CCNTPCL 

102 

5S 
1% 

515 
229 

51 

497 
254 

59 

502 

255 

46 

PSSF 
VS 
N    (0,1) 

10% 
1? 

3586 
3005 
1893 

4057 
3591 
2495 

4542 
4256 
3453 

•CCNTRCL'     IS    THE    TEST    CF    UNIFORMS, 1)    VS    UNIFORM(0,1) 

•RSSF'     IS    THE    RANDOM    STABILIZED    STANDARD    FORM 
*LPH*    =    1.00,     EETA    =    C.75 


NUMeER    OF    HYPOTHESIS    REJECTIONS     IN    5000    REPLICATIONS 


60 


CHI     SQUARED    TEST 

20    PTS                  2C    PTS  50    PTS 

10?                         368                         447  505 

CCNTROL          5%                         2CC                         241  222 

1?                            43                            58  39 


368 
2CC 

43 

19  41 

1495 

704 

PSSP  10*  1941  2S1 1  4C97 

VS  5?  1495  2429  361C 

I**    (  CU)  1*  704  1535  2762 


KOLMOG-SMIR    TEST 

20    PTS  3C    PTS  50    PTS 

10*          475  504  £24 

CCNTPCL    5?          242  237  264 

1*           47  45  47 

RSSF            10*                      2017  264C  3726 

VS                  5*                       1331  1930  3089 

N    (C»l)          1*                         546  904  1762 


4  75 

242 

47 

50  4 

237 

45 

20  17 

1331 

546 

264C 

1930 

904 

FCUTZ    FN 

TEST 

2  0    PTS 

30    DTS 

461 
226 

45 

513 

248 

54 

2895 
2227 
12C9 

3212 
2712 
1595 

50    PTS 

10*  461  513  490 

CCNTPOL         5*  226  248  246 
1* 

RSSF            1C*  2895  3212  4C51 

VS                 5*  2227  2712  3542 

N    (Oil)          1*  12C9  1595  2426 


•CCNTPCL*     IS    THE    TEST    CF    UNI  FORM ( 0 » 1 )    VS    UNIFCPN(Cl) 

•PSSF'    IS    1HE    RANDOM    STAEILIZEC    STANDARD    FORM 
ALPFA    =    1.20,    eETA    =    C.75 


M'NBEfl    CF    HYFCTHESIS    REJECTIONS     IN    50CO    REPLICATIONS 


61 


CHI    SQUARED    TEST 

20    PTS  3C    PTS  50    PTS 


CONTROL 

53! 
1? 

4C0 

2  43 

48 

407 

225 

51 

501 

256 

46 

RSSF 
VS 
N    (0,1) 

10? 

5? 
It 

1628 

1240 

509 

2451 
1965 
U26 

3633 
3C99 
211C 

CCNTFOL 

10? 

fZ<9 

li 

FSSF 
VS 

1*     (0»U 

10? 
5? 
1? 

KOLMOG-SNT*    T£ST 

20    PTS  3C    PTS  50    PTS 

476 

243 

50 


3022 
2330 
1145 


50     FTS 

1C  49C 

CONTROL         5?  266  267  242 

1?  52  60  56 

RSSF            10?  2426  2722  3367 

VS                 51  1753  2067  2727 

N    (0,1)         1%  825  1C58  1614 


491 

111 

48 

509 
244 

54 

1622 

10  50 
379 

2104 

1465 

588 

FOUTZ    FN 

TEST 

20    PTS 

3G    PTS 

522 

266 
52 

490 

267 

60 

2426 
1753 

825 

2722 

20  67 
1C58 

•CCNTROL1     IS    71-E    TEST    OF    UNIFQRM(0,1)     VS    UNI  FORM  (0,  1 ) 

•RSSF*     IS    THE    RANDOM    STAETLIZEC   STANDARD    FCPM 
ALPHA    =    1.60,    3ETA    =    0.75 


K'NBEF    OF    HYPOTHESIS    REJECTIONS    IN    5000    REPLICATIONS 


62 


CCNTPCL 

10* 
5? 
1* 

FSSF 
VS 
N  (0,1) 

10* 
5* 
1* 

CHI     SCU/!RED  TEST 

2C    PTS  30    PTS  50    PTS 

405  4C7  469 

?38  229  215 

54  49  46 


1366  1989  3088 

987  1490  2487 

341  704  I4ic 


KOLN0G-SMIR    TEST 

20    PTS  3C    PTS  50    PTS 


CONTROL 

10* 
5? 
1* 

501 

246 

57 

505 

268 

53 

483 

245 

46 

RSSF 

VS 

N  (0,1) 

10* 
5* 
1* 

1273 
797 

234 

FOUTZ  FN 

1558 
936 

308 

TEST 

2246 
1465 

508 

20  PTS 

30  PTS 

50  FT 

CONTROL 

1C* 

5* 
1* 

497 
242 

58 

457 

227 

48 

478 
242 

48 

RSSF 
VS 
h    (CD 

10* 
5* 
1* 

2102 

1479 

647 

2256 

16  15 

76  C 

2745 
2CE1 

1076 

•CCNTROL*     IS    THE    TEST    OF    UNIFGR^(0,1)     VS    UNIF0RM(0,1) 

•RSSF'     IS    THE    RANDOM    ST/*eiLTZEC    STANDARD    FCRM 
ALPHA    =    1.9C,     EETA    =    0.75 


f^LMBEf:-   OF    HYPCTHESIS    REJECTICNS    IN    5000    REPLICATICNS 


63 


CCN1FQL 


1C2 
5% 
1% 


CHI     SQUARED    TEST 
20    PTS  3C    PTS 


395 

213 

35 


441 

238 
53 


50    PTS 

495 

243 

62 


PSSF 
VS 
h    (0,1) 


10$ 

IS 


1793 

13<95 
554 


292C 
2386 
1413 


4109 
3658 
2707 


KGLMGG-SMIR    TEST 
20    PTS  30    P^S 


50    FTS 


CCNTFCl 

io-« 

5% 
1% 

493 

241 

50 

504 

231 

43 

4S7 

238 

43 

PSSF 
VS 
N  (CD 

1C1 
1% 

1640 

10  04 

295 

FOUTZ  FN 

2203 

139  2 

421 

TEST 

3325 

2434 

939 

20  PTS 

30  PTS 

50  PT 

CCNTPCL 

10? 

526 

260 
54 

488 
247 

58 

48  1 
223 

50 

PSSF 
VS 
N  (Ctl) 

1C3 
5% 
1% 

3530 
29  33 
1849 

4001 
3501 
2405 

4577 
4284 
34C5 

•CCNTftCL'     IS    THE    TEST    OF    UNIFORM(0,1)    VS    UNIFORMtOtlJ 

•FSSF«     IS    THE    RANDOM    STABILIZED    STANCARD    FORM 
JLPHA    =   1.00  t    BETA    =    1.00 


MMeEP    OF    HVPCTHESIS    REJECTIONS    IN    5000    REPLICATICNS 


64 


CHI     SQUARED    TEST 


20    PTS 

3C    PTS 

50    PTS 

CONTROL 

5% 
1? 

363 

218 

45 

428 

226 

44 

518 

243 
55 

RSSF 
VS 
N    <C,1  ) 

1C1 

5% 
1% 

2651 
2291 
1344 

2826 
3457 
2572 

4682 
4472 
3935 

KOLMOG-SMIR     TEST 


20    PTS 

30    PTS 

50    PTS 

CCNTPCL 

loss 

5% 

1% 

526 

243 

51 

533 

263 

48 

524 
262 

70 

PSSF 
VS 
N    <C,1) 

10* 
5% 
1% 

2937 
23  26 
1245 

3776 
3235 

2035 

4636 
4378 
3452 

FCUTZ    FN   TEST 


20    PTS 

30    PTS 

50    PTS 

CONTROL 

10* 

52 
1* 

488 

259 

54 

477 

216 

53 

500 

236 

44 

PSSF 

VS 

N     (0,1) 

102 

5% 
1% 

2  9  62 
2260 
1237 

3429 
2313 
1705 

4162 

3660 
2550 

•CONTROL*     IS    THE    TEST    OF    UNIF0RM(0,1)     VS    UNIF0RM(0,1) 

•PSSF»     IS    IhE    RANDOM    STABTLIZEC    STANDARD    FOP^ 
AIPH£    =    1.  20,     EETA    =    1,00 


NUMBER    OF    KVPC7HE5I5    REJECTIONS     IN    5000    REPLICATIONS 


65 


CHI     SQU4R  50    TEST 

20    PTS  3C    PTS  50    FTS 

10?                         412  443  472 

CCNTFOL         5%                        220  237  220 

1%                            30  50  46 

RSSF            1C3                      1868  2768  3S56 

VS                 5%                      1487  2294  3541 

N    (C,l  )         1?                        701  1427  2578 


KCLMQG-SMIR     TEST 

20    PTS  3C    PTS  50    PTS 


CCNTPCL 

10? 
5? 
1? 

513 

268 

56 

477 

251 

56 

514 

251 

49 

RSSF 
VS 
N     (0,1) 

10? 

5% 

1968 

1370 

586 

POUTZ    Ffs 

2603 

1903 

905 

TEST 

3652 
30  16 
1765 

20    PTS 

3C    PTS 

50    FT 

CGNTPGL 

10? 
53 
1? 

507 

25C 
58 

504 

240 
45 

5C5 

252 

58 

RSSF 
VS 
N     (C.I) 

10? 

5S 
1? 

2414 

1786 
853 

2792 
2120 
1123 

3409 
2771 
1625 

•CCNTRaL'     IS    THE    TEST    CF    UMIFORM(0TD    VS    UNIFORM  (Ct  1  ) 

•RSSF*     IS    THE    RANCOM    STABILIZED    STANDARD    FQPM 
jiLPHJS    =    1.60,     EETA    =    l.OC 


M.MEEP    Of    HYPCTHESIS    REJECTIONS     IN    5000    P.  EPL ICAT ICNS 
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CHI     SQUARED    TEST 

20    PTS  3C    PTS  50    P7S 

1C*                        374  431  525 

CCNTPOL         51                       231  232  263 

1%                           48  57  58 

FSSF     10*         1339  1S53  31C7 

VS       5%                          977  1438  2466 

N  (0,1)    1*          347  708  1459 


KOLMOG-SMIR    TEST 

20    PTS  3C    PTS  50    PTS 

IZ%                        494  481  516 

CCN7PCL    5%                        229  240  241 

1%                           54  47  43 

PSSF            10?                      1322  1625  23CE 

VS                  53                         8C8  978  1533 

N    (CI)         1$                        235  302  540 


FOUTZ    FN   T  EST 

20    PTS  30    PTS  5C    PTS 


CCNTFCL 

10* 
5% 
1* 

525 

236 
50 

523 

271 

67 

520 

296 

5* 

PSSf 
VS 
N    (CD 

1CS 

5% 
1% 

2086 

1409 

6  20 

2280 

1645 

715 

2758 

2075 
1078 

•CCNTFCL1     IS    THE    TEST    OF    UN  IFQRMC  0  ,1 )    VS    UNXFCRM  (0  1 1) 

•FSSF»    IS    THE    FANCCM    STAEILIZEC    STANDARD    FORI* 
ALPHA    =    1  .90  ,    EETA    -    1.00 


M'MeEP    CF    HVPCTHESIS    REJECTIONS     IN    50C0    REPLICA1ICNS 
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CHI     SQl/RED    TEST 

2C    PTS  3C    PTS  50    PTS 


CCNTFCL 

10? 
5* 
131 

388 

2  29 

39 

469 

238 

48 

4E6 
224 

43 

MXNCRN 

\S 
N     (C»l) 

1C? 
5? 
1% 

12  72 

909 
329 

1925 

1407 

664 

30C4 
2361 

1337 

KOLMOG-SMIR     TEST 

20    PTS  30    PTS  50    PTS 


CCNTRCl 

10? 
5? 
1? 

529 

254 
42 

513 

263 

46 

5  10 

249 

64 

MIXNCRM 

VS 
N    (Ctl) 

10? 
5? 
IS 

1235 
7  13 
208 

FGUTZ    FN 

1537 
931 
28  5 

TEST 

2133 

1369 
449 

20    PTS 

3C    PTS 

50    PT 

CCNTPOL 

10? 
15 

508 

245 

£0 

476 

222 

38 

523 

255 

6C 

MXNCPM 

VS 
N     (0,11 

10? 
5? 
1? 

1982 

1336 

533 

2191 

1576 

672 

26  10 

1995 

977 

•CCNTRai1     IS    THE    TEST    CF    UNIF0RM(0,1)     VS    UNIFORMCl) 

•MXNCRN*     IS    ThE    MIXED    NORMAL    DISTRIBUTION 
C-  C   *N(.Otl)    ♦    1.0C*N<  C,2  ) 


K'NBEF    GF    HYPOTHESIS    REJECTIONS    IN    5000    REPLICATICNS 
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CCMFCL 


10? 
5? 
1? 


CHI    SQUARED    TEST 
20    PTS  3C    PTS 


4C5 

233 

29 


411 

205 

34 


5C    PTS 

498 

234 

46 


NIXNCRM       10? 

VS  5? 

N    (0,1)         12 


2513 

20  88 
1028 


3558 
3C8  0 
2052 


462C 
4342 
36C9 


CCNTSOl 


CCNTPQL 


10? 
5? 
1? 


MXNCPM       10? 

VS  52 

N    (0,1)         1? 


10? 
5? 
It 


MXNCFM       10? 

VS  5? 

h    (0,1)  1? 


KOLMOG-SMIR    TEST 
20    PTS  3C    PTS 


475 

238 

46 

463 

229 

35 

2188 

1420 

512 

2827 

1974 

727 

FOUTZ  F  h 

TEST 

20  PTS 

3C  PTS 

522 

236 

47 

454 

217 

41 

3338 
2724 
1657 

3712 
3174 

2064 

50  PTS 

50C 
242 

38 


39  6  5 
3161 
1515 


50  PTS 

481 

226 

52 


4346 

3935 
2970 


•CONTROL*     IS    THE    TEST    CF    UN  I  FORM  (0  ,  1)    VS    UNIFORM(Cl) 

•MXNCPW     IS    THE    MIXED    NORMAL    DISTRIBUTION 
C.C    *M.0,1)    +    1.00*N(0,3) 


f>UNBEF    OF    HYPOTHESIS    REJECTIONS     IN    5000  R  E?  L  IC  A  TIQNS 
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CONTROL 


1GS 
5% 
1% 


CHI     SQUARED    TEST 
20    PTS  3C    PTS 


407 

2  34 

48 


425 
232 

46 


50    PTS 

511 

252 
5C 


MIXNQFM       102 
VS  5% 

r>   (0,1)       i? 


3427 
3006 
17£1 


4423 
4154 
3377 


4935 
4863 
46  39 


CCNTPQL 


CCNTFOL 


1C3 
535 

12 


IM>NCRN       1C* 

VS  5% 

N    (Ctl)         1? 


10? 
5% 
1% 


MIXNCPM       IZ% 

VS  5* 

N    (C,l  )         1? 


KOLMOG-SMIR    TEST 
20    PTS  3C    PTS 


510 

261 

59 

490 

241 

48 

2976 

2113 

890 

3799 
30  25 
1471 

FOUTZ  FN 

TEST 

20  PTS 

30  PTS 

518 

242 

68 

501 

268 

66 

4131 
3686 
2653 

4502 
4195 

3380 

50    PTS 

511 

250 

47 


4717 
42  66 
2792 


5C    PTS 

528 

276 

62 


4E65 
4742 
4272 


•CCNTPCL*     IS    THE    TEST   OF    UNIF0RM(C,1)    VS    UNI  FCP-V  (C ,  1> 

•VIXNCRN*     IS    THE    MIXED    NORMAL     DISTRIBUTION 
0.0    *N<  .Otl  )    +    1.00*N(0t4l 


MM8EP    OF    FVPCTHESIS    REJECTIONS     IN    5000    REPLICATIONS 
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CHI  SQUARED  TEST 

20  PTS  30  PTS  50  PTS 

10?                        382  422  510 

CCNTFOL         53                         215  137  254 

i%                           41  38  46 

M*KR*      10?                        472  542  713 

VS                  5?                         2S3  308  370 

I*    (Ctl)          1?                            29  66  74 


KOLMOG-SMIR    TEST 

20    PTS  30    PTS  50    PTS 

1CS          488  537  525 

CCNTFCL    5%                        238  250  265 

1?           46  50  53 

MIXNCRP   1C?          587  667  695 

VS       5%                        287  328  367 

N  (CD    U           54  70  84 


FOUTZ  FN  TEST 

20  PTS  30  PTS  50  PTS 

10%                        468  480  517 

CCNTPCt         5%                         220  258  247 

1*                            46  64  56 

MIXNCRN      10?                         799  796  861 

VS                 5%                         429  426  477 

N    (0»1)          II                        105  107  136 


■CCNTRCL1     IS    THE    TEST    CF    UNIF0RM<0,1)    VS    UNIFORM(Ctl) 

•MIXNQRM*     IS    THE    MIXED    NCPMAL    DISTRIBUTION 
C.7C*N(.0»1)     +    C.30*N<0,2) 


MMEEP    OF    HYPOTHESIS    REJECTIONS     IN    5000    REPLICATIONS 


71 


CCNTFOL 


10* 
5% 
1* 


CHI     SQU4RE0    TEST 
20    PTS  3C    PTS 


421 

245 

49 


457 

244 

52 


5C    PTS 

510 

221 

36 


NIXNCRM       1C* 

VS  5* 

N    (CD         1% 


551 

345 

82 


7  27 
438 
131 


1C4  7 

602 
185 


CCNTPOL 


CONTROL 


10* 
5* 
1* 


MIXNORM       10? 

VS  5? 

N    (0,1)         I* 


10* 
5* 
1% 


MIXNCRM       10* 

VS  5? 

h    (0,1)  1* 


KOLMCG-SMIR    TEST 
20    PTS  30    PTS 


490 

251 

46 

504 

256 

58 

6  46 

337 

81 

745 

405 
93 

FOUTZ    F!^ 

TEST 

20    PTS 

3C    PTS 

5  19 
24C 

55 

494 

242 

46 

1029 
591 
172 

1139 
671 
193 

50    PTS 

527 

23< 

42 


845 

447 

€3 


50    PTS 

507 
256 

51 


1249 
77£ 

262 


•CCNTRCL*      IS    THE    TEST    CF    UNIFORMS, 1)    VS    U  NI  FORN  (  C,  1  > 

•MlXNORf"     IS    THE    MIXED    NORMAL    DISTRIBUTION 
0.7C*N(.0,1)    +    0.30*N<0,3) 


MJV8EF    CF    HYPCTHESIS    REJECTIONS     IN    5000   REPLICATIONS 
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CCNTPGL 


1C3 
5% 
13 


CHI     SQUARED    TEST 
20    PTS  3C    PTS 


400 

237 

45 


435 
223 

50 


50    PTS 

515 

240 

52 


MXNCRV 

vs 
r>   (0,1 ) 


10? 

13 


626 
4C5 
109 


S10 
564 
175 


1320 
336 


CCNTFCL 


CCNTPCL 


1C3 
5% 
1% 


fIXNCRN       1C? 

VS  5% 

N    (C,l)         1? 


10? 
5% 
IS 


NIXNCRN       10? 
\j  c  52 

N    (6,1)  1? 


KOLMOG-SMIR     TEST 
2C    PTS  30    PTS 


508 

259 

54 

531 

259 
51 

747 
407 
115 

853 
441 
122 

FOUTZ  FN 

TEST 

20  PTS 

30  PTS 

492 

229 

61 

478 

222 

47 

1281 
774 
238 

1418 
888 
222 

50  PTS 

5CE 

248 

54 


10C1 
532 
148 


50  PTS 

511 

241 

42 


1666 

1121 

4C1 


•CCNTRCL»  IS  THE  TEST  CF  UNIF0RM(0,1)  VS  UNIFORM(Otl) 

•MIXNORM*  IS  THE  MIXED  NCRMAL  DISTRIBUTION 
C.7C*N(.0,1)  +  0.30*N<0,4) 


M.MEEP  GF  HYPCTFESIS  REJECTIONS  IN  5000  REPLICATICNS 
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CHI    SQUARED    TEST 
20    PTS  30    PTS 


50    FTS 


CCNTfCL 

10? 
5? 
1? 

380 

221 

40 

431 

245 
59 

5C7 

220 

56 

MXNCRM 

VS 
N    (0,1) 

ic? 

IS 

388 

214 

40 

KOLMOG-S^n 

484 
274 

63 

TEST 

545 

275 

62 

20    PTS 

30    °TS 

50    PT 

CCNTPCL 

10? 

5? 
1* 

468 

220 
53 

480 

236 

46 

481 
232 

28 

NIXNCRM 

VS 
N    (Otl) 

10? 

5? 
1? 

519 

258 

58 

539 

277 

55 

556 

286 
53 

CCNTPOL 


10? 
5% 
1? 


FOUTZ    FH   TEST 

20    PTS  3C    °TS 


487 

234 

42 


498 

242 

49 


50    PTS 

502 

265 

54 


*IXNCFM       10? 

VS  5? 

N     (0»1)  1? 


684 

352 

63 


689 

354 
92 


687 

400 

96 


•CCNTRGL'      IS    THE    TEST    CF    UNIFQRM(0,1)    VS    UNIFORM  (C,l) 

•*!XNCPM»     IS    THE    MIXED    NORMAL    DISTRIBUTION 
C.€C*M.0»1)    ♦    0.20*N<0,2) 


t^UNBEF    OF    HYPOTHESIS    REJECTIONS    IN    5000    REPLICATICNS 
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CfNTFOL 


1CS 

1% 


CHI    SQUARED   TEST 
20    PTS  20    PTS 


417 

218 

43 


405 
211 

38 


5C    PTS 

492 

243 

42 


MXNCRM       1C1 

VS  5% 

IS    (C,l  )         1* 


477 

281 

48 


544 

306 

64 


736 

364 

98 


CCNTRCl 


CONTROL 


10? 
IS 


MIXNCRM       1C9? 

VS  5? 

N-(C,1)         1% 


10? 
5? 

1? 


MXNCPM       10% 

VS  5% 

N     (Ctl)  1% 


KOLMOG-SMIR    T FST 
20    FTS  30    PTS 


5CC 

267 

47 

469 

222 

46 

587 

3  17 

69 

590 

297 

59 

FOJTZ    FN 

TEST 

20    PTS 

30    PTS 

495 

234 

49 

46? 

213 

44 

839 
475 

118 

819 
454 
130 

50    PTS 

527 

26C 

46 


676 

342 

16 


50     PTS 

479 
245 

48 


945 
522 
142 


•CCNTRCL'     IS    THE    TEST    CF    UNIFORM(Otl)     VS    UNIFORMED 

•MIXNCRV     IS    THE    MIXED    NORMAL     C ISTR I8UTI0N 
C.£C*N<.C,1)     +    O.20*N<0,2) 


NUMBER    OF    HYPOTHESIS    REJECTIONS    IN    5000    REPLICATIONS 
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CHI     SQUARED    TEST 

20    PTS  3C    PTS  50     PTS 

10?          402  450  539 

CCMPCt    5?          225  242  247 

1?                            43  49  6C 

VIXNCRM       1C?                         507  637  9C2 

\S                 5%                         310  361  496 

N    (CI)          II                            74  89  138 


KOLMOG-SMR     TEST 

20    PTS  30    PTS  50    PTS 


CCNTPCL 

10? 
5% 
1? 

511 

253 

51 

502 

2  52 

43 

542 

270 

£3 

MXNCRM 

VS 
N    (0,1) 

10? 
5? 
1? 

622 
3  29 

85 

FCUTZ    FN 

661 

360 
69 

TEST 

713 

416 
S3 

20    PTS 

30    PTS 

50    PT 

CCNTPCL 

10? 
5? 
1? 

485 

238 

48 

496 

243 

48 

514 

259 

c  c 

MIXMOFM 

VS 
M    (0,1) 

10? 
5% 
1? 

958 
548 
171 

1011 
583 
174 

1168 
721 

24* 

•CONTROL'     IS    THE    TEST    CF    UNIFCRM(OU)    VS    UNIFQRMO,!) 

•MI>NCRN«     IS    THE    MIXED    NCRM4L    DISTRIBUTION 
C.£C*N(.0fl)    +    0.20*N(C,4) 


KNBEF    CF    HYPCTHESIS    REJECTIONS    IN    5000    REPLICATICNS 
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CONTROL 


OS 

1% 


CHI    SQUARED    TEST 
20    PTS  3C    PTS 


420 

231 

40 


470 

240 

51 


50    PTS 

47C 

23£ 
44 


M  >NGPM 

VS 
1^     (0,1) 


5% 
1% 


41C 

228 

35 


478 

241 

48 


525 

247 

46 


CCNTPOL 


CONTROL 


1C? 
53 

13? 


M>NCRM       101 

VS  5? 

h    (CD         13 


10% 
52? 
13 


MXNCRf-      ICS 
N    (CD        13 


KOLMOG-SMR    TFST 
20    PTS  3C    PTS 


499 

244 

54 

490 

245 

44 

531 

264 
57 

512 

257 
47 

FOUTZ  FN 

TES 

T 

20  PTS 

30  PTS 

493 
245 

43 

490 

256 

62 

577 

282 

66 

' 

583 

310 

73 

50    PTS 

519 
251 

55 


535 

266 

64 


50    PTS 

4EC 
229 

47 


587 

2<C 

77 


•CCNTPCL*     IS    THE    TEST    OF    UNI  FOR  M<0  ,1 )    VS    JNlFORMOtl) 

•MIXNORM*      IS    Tl-E    MIXED    NORMAL    CISTRIBUTICN 
0.9Q*NC0tl)    +    0.10*N(0,2) 


NUNeEP    OF    I-YPCTHESIS    REJECTIONS    IN    5000    REPLICATICNS 
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CHI    SQU4R  ED    TEST 

20    PTS  30    PTS  50    FTS 

10?                        432  416  466 

CCNTRCl         5%                          246  214  215 

1?                           43  41  4C 

MXNCRI*       10?                         456  431  537 

VS                 53                         262  233  2  58 

N    (CtlJ         1?                           51  43  53 


K0LM0G-SMIR    TEST 

20    PTS  30   PTS  50    PTS 

10*          5C9  495  5CS 

CCNTfOL    51          271  239  252 

1%                            55  46  55 

VIXNCPM       1C3                         553  526  557 

VS                 5%                          302  261  274 

N    (0,1)          I?                            67  53  64 


5C9 

271 

55 

495 

239 

46 

5  53 

302 

67 

526 

26  1 

53 

FQUTZ  FN 

TEST 

2  0  PTS 

30  DTS 

516 
271 

58 

471 

239 

51 

66C 

353 
95 

653 

357 

88 

50    PTS 

10  a  520 

CCNTPQL         5?  271  239  252 

1%  58  51  47 

PIXNCRM       10?  66C  653  675 

VS                 5?  353  357  356 

N    COtl)         1?  95  88  ES 


•CCNTRCL1     IS    THE    TEST    OF    UNIF0RM(0,1)    VS    UNIFCRM(Otl) 

•M>NCPVr     IS    THE    MIXEC    NORMAL    DISTRIBUTION 
C.5C*N(.0,D     ♦    0.1G*N(0»3) 


NUMEER    OF    HYPOTHESIS    REwECTICNS    IN    5000    REP  LI  CAT  ICf^S 


78 


CONTROL 


IC35 

1% 
1% 


CHI    SQUARED    TEST 
20    PTS  30    PTS 


3S7 

240 
31 


439 

244 

53 


50    PTS 

509 

246 

55 


MXNCRN       1C2 

VS  5% 

N    (CD         1% 


419 
260 

36 


508 

266 

65 


551 

3C2 

74 


KOLMQG-SMIR    test 
2C    PTS  30    PTS 


5C    PTS 


CCNTfOL 

1035 
5% 

155 

4  84 

240 

56 

516 

247 

47 

466 

239 

49 

HIXNCRN 
VS 

N    <C,1  ) 

1C? 
5% 
1? 

525 

268 
64 

FOUTZ    FN 

TEST 

584 

282 

61 

54€ 
217 

6C 

2C    PTS 

30    PTS 

50    PT 

CCNTPCL 

10? 
53! 
1% 

514 
?47 

48 

511 

253 

50 

49  1 
245 

58 

MIXNCRM 

VS 
IS    <0,1) 

10? 
5? 
1? 

736 

385 

81 

707 

384 
105 

747 
422 
112 

•CCNTRCL*     IS    THE    TEST    CF    LNIFQRM<0,1)    VS    UNIFCRMCCtl) 

•HXNCRN1    IS    THE    MIXEC    NCRMAL    CISTRIBUTION 
C.90*NM,0,1)     +    O.10*N(0,4) 


NUMBER    OF    HYPOTHESIS    REJECTICNS    IN    5000    REPLICATICNS 
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CHI     SQU/5R  ED    TEST 

20    PTS  30    PTS  50    FTS 

10%          386  417  4*4 

CCNTPCL    5%                         225  241  238 

1%                             39  44  46 

HXNCRP       103                       1459  2008  3146 

VS                 5%                      1093  1533  2556 

is    (Ctll         1%                         407  793  1568 


KCLMOG-SMIR    TEST 

20    FTS  3C    PTS  5C    PTS 


CCNTFOL 

10% 

5% 
1? 

493 

267 

52 

48  3 

246 

53 

487 

238 

36 

MIXNCRN 

VS 
N  (Ctll 

1C2 
5? 
1? 

2179 

1606 

714 

FOUTZ  FN 

2806 
2181 
1148 

TEST 

3756 
3226 

2043 

2C  PTS 

30  PTS 

50  PT 

CCNTPCL 

10? 
5* 
1% 

535 

250 

54 

495 

246 

43 

499 

248 

41 

MIXNGRM 

VS 
N  (0,1) 

10? 
5% 
1% 

18C9 

1198 

509 

1994 

1428 

602 

2492 

1825 

865 

•CCNTRCL'     IS    THE    TEST    OF    UNIFORM(0,U    VS    UNIFORM  (C,1J 

«MI>NOPM«     IS    THE    MIXED    NCRMAL    DISTRIBUTION 
C.7C*N(.5»1)     ♦0.3C*N(C,3) 


NUMBER    OF    HYPOTHESIS    REJECTIONS    IN    5000    REP  LICATI  CI^S 
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CHI    SQUiR  ED    TEST 

20    PTS  30    PTS  50    FTS 

10%                         377  449  512 

CCNTRCL         5?                         205  239  257 

IS                             46  46  44 

NIXNCRM      10?                      1557  2112  3304 

VS                 53?                      1170  1563  2660 

N    (0,1)          1?                        473  750  164C 


KOL^CG-SMIR    TEST 

20    PTS  30    PTS  50    PTS 

10?                         482  481  538 

CCNTRCL         5%                        229  260  271 

1?                            45  69  56 

PIXNCRM       10?                      245C  3101  4056 

VS                 5?                      1855  2510  3596 

N    (0,1)          1?                         9C1  1344  2400 


482 

229 

45 

481 

260 

69 

245C 

1855 

9C1 

3101 
2510 
1344 

FOUTZ    FN 

TEST 

20    PTS 

3C    PTS 

479 

233 

50 

512 

262 

53 

1781 

12C9 
462 

1917 

1369 
573 

50    PTS 

1C?  479  512  531 

CCNTRQL         5?  223  262  255 

1?  50  53  53 

HXNCRN      10?  1781  1917  2405 

VS                 5?  12C9  1369  1774 

N    (0,1)         1?  462  573  842 


•CCNTRCL'     !S    THE    TEST    CF    LM  FORM  <0, 1)    VS    UNIFCRM(Cl) 

«*IXNCRN«     I<    THE    MIXEC    NCRMAl    CISTRIBUTICN 
0.€C*M.5tl)    ♦    0.20*N(0,3) 


MJN8EP    CF    HYPOTHESIS    REJECTIONS     IN    5000   REPLICATICNS 
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CHI    SQLUREO    TEST 

20    PTS  3C    PTS  50    F7S 

103?                         389  439  490 

CCNTFCL         5*                       221  251  221 

1%                          38  55  49 

MXNCRM   1C?         1568  2228  3377 

\S       5%                      1177  1736  2832 

N  (CD    1*          507  951  1796 


KOLMOG-SMR     TEST 


20    PTS 

30    PTS 

50    PTS 

CCNTRCL 

10% 
5% 
1% 

535 

258 

55 

489 

253 

54 

490 
24c 

35 

VIXNCRM 

VS 
N    (0,1) 

10* 
5% 
1% 

2467 
1865 

857 

3194 
2593 
1438 

41C5 
3659 
2511 

FCUTZ    FN   TEST 


20    PTS 

30    PTS 

50    PTS 

CONTROL 

10% 

5% 
1% 

502 

2  52 

62 

496 

250 

51 

495 
23c 

c  -z 

MXNCFP* 

VS 
N     (0,1) 

10% 
5% 
1% 

183  8 

1220 

486 

2156 

1551 

701 

2665 

1978 

979 

•CONTROL*     IS    THE    TEST    CF    UNIF0RM(0,1)     VS    UNIFQRMOtl) 

■NDNCRN*     IS    THE    MIXED    NORMAL    DISTRIBUTION 
C.€C*NU5tl)    ♦    0.20*N(C,4) 


K'NBEF    CF    HYPCTHESIS    REJECTIONS    IN    5000    REPLICATIONS 
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CCNTFCL 


10? 
52 
12 


CHI  SQUARED  TEST 
20  PTS       3C  °TS 


431 

240 

40 


413 

200 
42 


50  PTS 

526 

249 

53 


FEAGSCN 

VS 
N  (C»l) 


10? 
5? 
1? 


465 

258 
46 


478 
242 

59 


tit 

302 

71 


CCNTFCL 


CCNTFCL 


102 
5% 
1? 


FEARSCN   10? 

VS       52 

N  (0,1)    1? 


1C? 
52 
1? 


FEAPSCN   102 

VS       52 

N  (CD    1? 


KOLM0G-SMIR  TEST 
20  PTS       3C  PTS 


528 

275 

45 

478 

242 

49 

623 
339 

66 

571 

310 

70 

FOUTZ  FN 

TEST 

20  PTS 

3C  PTS 

512 

248 
56 

458 

207 

48 

481 
241 

60 

428 

196 

36 

50  PTS 

525 

283 

54 


673 

371 

95 


50  FT? 

53C 

Itt 

58 


4<;7 

234 

57 


•CCNTFCL1  IS  THE  TEST  OF  UNIP0RM(0,1)  VS  UNIFORM (0 ♦ 1) 

•PEARSCN*  IS  THE  PFAPSON  DISTRIBUTION 
8ET£  1  =  0.0  ,  BETA  2  =  2.30 


NLMBEF  OF  HYPCTHESIS  REJECTIONS  IN  5000  REPLICATICNS 


83 


CHI     SQUARED    TEST 


20  PTS 

3C  PTS 

50  PTS 

CCNTFCL 

10? 
5? 
1? 

435 

227 

49 

443 

241 

52 

541 

247 

49 

FEAPSCN 

VS 
N  (C,l) 

10? 
5? 
1? 

432 

243 

43 

445 
247 

55 

532 

241 

52 

CCNTFCL 

10? 
5? 
1? 

FEAPSCN 

VS 
N  (0,1) 

10? 
5? 
1? 

CCNTFCL 


10? 
5? 
1? 


FEAPSCN   10? 
VS       5? 

h  ten       i? 


KOLNCG-SNIR    TEST 
20    PTS  3C    PTS 


489 
238 

55 

528 

252 

50 

496 

242 

61 

543 

274 

54 

FOUTZ  FN 

TEST 

20  PTS 

3C  PTS 

536 

275 

61 

467 

227 

34 

523 
269 

58 

441 

217 

33 

50  PTS 

489 

245 

64 


502 

262 

tl 


50  FTS 

487 

240 

43 


467 

228 

43 


•CCNTFCL'     IS    THE    TEST    OF    UMFCRM(0»1)     VS    UNIFORM(Cl) 

•PEARSON'     IS    THE    PEAPSON    DISTRIBUTION 
BETA    1    =    0.0    ,     BETA    2    =    2.80 


NLMBEP    CF    HYPCTHESIS    REJECTIONS     IN    50C0    REPLICATIONS 
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CCNTFCL 


1C* 

5? 
1? 


CHI  SQUARED  TEST 
20  FTS       3C  PTS 


4C8 
230 

35 


390 

214 

50 


50  PTS 

452 
213 

52 


FEAPSCN 

VS 
N  (CI) 


1C? 
5? 
1? 


756 
495 
121 


992 
634 

185 


1576 
973 
343 


KCLVCG-SMIR  TEST 


20  PTS 

3C  PTS 

50  FTS 

CCNTFCL 

10? 
5? 
1? 

492 

252 

55 

515 

255 

49 

520 

264 

52 

FEAPSCN 

VS 
N  <C,1) 

10? 
5? 
1? 

942 
5  46 

158 

1161 
693 

192 

1498 
944 
294 

FOUTZ  FN  TEST 


20  PTS 

3C  PTS 

50  PTS 

CCNTFCL 

10? 
5? 
1? 

5C5 

2  23 

50 

478 

261 

56 

5C4 
255 

48 

FEAPSCN 

VS 
1^  (CD 

10? 
5? 
1? 

853 
469 
129 

1400 
884 
311 

2487 

1785 

799 

•CCNTFCL*  IS  THE  TEST  CF  UNIFCRM(0,1)  VS  UNIFORMS, 1) 

•PEARSCN*  IS  THE  PEARSON  DISTRIBUTION 
BETA  1  =  0.01,  BETA  2  =  1.75 


NLfeEF  CF  HYPCTHESIS  REJECTIONS  IN  5000  REPLICATICNS 
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CHI  SQUARED  TEST 


20  PTS 

3C  °TS 

50  PTS 

CCNTFGl 

10? 
5? 

11! 

428 

239 

40 

410 

205 

43 

492 
52 

FEAPSCN 

VS 
H  (CI) 

10? 
53 
1? 

614 
375 

93 

744 
43  5 
120 

1082 
610 
193 

CCNTFGL 


10? 

5% 
13 


FEAPSCN   10? 

VS       5? 

N  (CD    13 


IC3 

CCNTFCl    5? 

1? 


PEARSON   10? 

VS       5? 

f>  (CD    1? 


KOLMOG-SVIR  TEST 
20  PTS       3C  PTS 


543 

290 

41 

470 
229 

50 

837 
482 
131 

857 
479 
131 

FOUTZ  FN 

TEST 

20  PTS 

3 

C  PTS 

511 

240 
55 

474 

214 

49 

558 

305 

82 

726 

405 
ICl 

50  PTS 

5C6 
274 

48 


1106 

648 
184 


50  PTS 

515 

259 

54 


1288 
774 
233 


•CCNTFCL1     IS    THE    TEST    OP    UNIFORM(Cl)    VS    UNIFCRM(Otl) 

•PEARSON*     IS    the    PEARSON    DISTRIBUTION 
BETA    1    =    O.OD    BETA    2    =    1.90 


M'MBEP    OF    HYPCTHESIS    REJECTIONS    IN    5000    REPLICATIONS 
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CHI    SQUARED    TP3T 
20    PTS  3C    PTS 


50    PTS 


CCNTFCL 

10? 
5? 
1% 

408 

23C 

48 

442 

237 

41 

510 

24! 

47 

FEAFSCN 

vs 

N  (C,l) 

10? 
5? 
1? 

458 

267 

51 

556 

321 

74 

702 

366 

96 

CCNTFCL 


CCNTFCL 


10? 
5? 
1% 


FEARSCN       1C? 

V c  5? 

N    (C,l)  1? 


10? 
5% 
1? 


PEARSCN       101 

\<S  5? 

N     (Ctl)  1? 


K0LM0G-SVIR    TEST 
20    PTS  3C    PTS 


49"? 

240 

55 

485 

272 

50 

581 

3C9 

74 

65C 

229 

73 

FOUTZ  FN 

TEST 

20  FTS 

30  PTS 

529 

241 

62 

479 

241 

56 

534 

269 

68 

554 

283 

66 

50  FTS 

5C2 

258 
65 


715 

3£9 

94 


50  PTS 


546 

285 

54 


666 

353 

81 


•CCNTFCL'  IS  THE  TEST  CF  UNIFORM(0,l)  VS  UNI  FORM  (0  , 1 ) 

•PEARSCN*  IS  THE  PEAFSON  DISTRIBUTION 
BETA  1  =  0.25,  BETA  2  =  3.20 


MMEEF  OF  HYPOTHESIS  REJECTIONS  IN  5000  REPLICATIONS 
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CCNTPOL 


10* 
5% 
1% 


CHI  SQUARED  TEST 
20  PTS       3C  PTS 


373 

22C 

55 


419 

230 

49 


50  PTS 

502 

211 

34 


FEAPSCN   10? 

V c       535 

N  (C,l)    1% 


739 
476 
118 


1C30 
664 
231 


1673 

1099 

431 


CCNTFGL 


CCNTFCL 


10? 
5? 
1% 


FEAPSCN   10? 

VS       5? 

h  <C,1)    1? 


1C3 
5? 
11 


FEARSCN   10? 

VS       5? 

N  <C,1)    1? 


KOLMOG-SMR  TEST 
20  PTS       3C  PTS 


516 

253 

48 

483 

227 

51 

790 
474 
121 

896 
528 
147 

FOUTZ  FN 

TEST 

20  PTS 

3C  PTS 

499 

237 

46 

481 

238 

49 

868 
471 
121 

1090 
660 
185 

50  FTS 

475 
231 

38 


1125 
735 
234 


50  PTS 

477 

253 

5C 


1516 
965 
344 


«CCNTFCl«     IS    1FE    TEST    OF    UNI  FOR  M(  C  ,  1 )  .  VS    UNIFORM(Cl) 

•PEARSCM     IS    THE    PEARSON    DISTRIBUTION 
EET*    1    =    0.50,    SETA    2    =    3.00 


M*BEF    OF    HYPC7HESIS    REJECTIONS    IN    5000    REPLICATICNS 
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CCNTFCL 


1CS 

5% 
1% 


CHI  SQUARED  TEST 
20  PTS       3C  PTS 


4C1 

236 

43 


428 

225 

39 


50  PTS 

22S 

54 


FEAPSCN   10? 
VS       5% 

r*  ecu      1? 


2440 

1960 

931 


3386 
2620 
1632 


4654 
432C 
3267 


KOLMOG-SMIR  TEST 


20  PTS 

3C  PTS 

50  PTS 

CCNTFCL 

1C3 

5? 
11 

482 

239 

47 

513 

254 
52 

5C9 

245 

61 

FEAPSCN 

VS 
h  (CD 

1C? 

5? 
1% 

1155 

730 
216 

1459 

920 
322 

20C1 

1428 

6C9 

CCNTFCL 


10% 
5? 

1% 


FOUTZ  FN  TEST 

20  PTS       30  PTS 


531 

250 

57 


492 

250 

51 


50  PTS 

55C 

277 

61 


FEAPSCN   1C? 

VS       5% 

N  (CD    1? 


2135 

1460 

58^ 


2837 
2156 
1049 


392C 
3309 
21C4 


•CCNTFCL1     IS    THE    TEST    CF    UNIFORM(0,1)    VS    UNIFORM<0»1) 

•PE/>RSCN»     IS    THE    PEAFSON    DISTRIBUTION 
BETA     1    =    l.OCt    BETA    2    =    3.40 


MMEEF    CF    HYPCTHESIS    REJECTIONS     IN    5000    REPLICATICNS 
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CHI  SQUARED  TEST 
2C  PTS       30  P^S 


50  PTS 


CCNTFCL 

10? 
5? 
1? 

411 

226 

34 

427 

238 

49 

475 

227 

34 

FEARSCN 

vs 
n  (CD 

10* 

5? 
1? 

1446 
1070 

4C3 

3151 
2546 
1368 

4513 
4039 
2895 

CCNTFCL 


10? 

5% 
1% 


KOLMQG-SMIR  TEST 
20  PTS       3C  PTS 


499 
249 

55 


495 

260 

64 


50  PTS 

5C4 

26  9 

6C 


PE4PSCN   1C? 

VS       5? 

N  (CD    1? 


1032 
649 
211 


1314 

832 
265 


1765 

1221 

493 


CCNTFCL 


10? 
5? 
1* 


FCUTZ  FN  TEST 

20  PTS       3C  PTS 


482 

251 

51 


481 

251 

45 


50  PTS 

485 

2C7 

43 


FEAFSCN   10? 
V<       5| 

N  (CI)    1? 


1578 

1009 

348 


2161 
1493 

633 


3119 
2354 
1254 


•CCNTFCL'  IS  ThE  TEST  CP  UNI  FORM  (  0  , 1 )  VS  UNIFORM(Cl) 

•PEARSCN*  IS  THE  PEARSON  DISTRIBUTION 
BET/!  1  =  1.00,  BETA  2  =  3.60 


NLVEEF  OF  HYPOTHESIS  REJECTIONS  IN  5000  REPLICATIONS 
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CCNTFCL 


1CS 
5? 
1? 


FEAPSCN   1C? 

VS       5% 

K  (C1J    1? 


CHI  SQUARED  TEST 
20  FTS       3C  PTS 


397 

23? 

41 


1086 
764 
255 


427 

233 
60 

2C09 

1463 

651 

50  PTS 

469 

234 

56 


3226 
2542 
136'' 


CCNTFCL 


CCNTFCL 


10? 
5? 
1? 


FEAFSCN   10? 

VS       5? 

N  <C,1)    1? 


10? 
5? 
1? 


FEAPSCN   10? 
V*       5? 

r  (c,d      i? 


KCLNCG-SMIR    TEST 
20    PTS  3C    PTS 


524 

267 

51 

501 

250 

56 

911 
546 
171 

1150 
732 
229 

FOUTZ  FN 

TEST 

20  PTS 

3C  PTS 

521 

?58 

50 

460 

248 

63 

1261 
764 
243 

1656 

1070 

387 

50  PTS 

499 

252 

62 


1526 

1022 

299 


50  PTS 

529 

260 

46 


2425 

17C6 

729 


•CCNTFCL1     IS    THE    TEST    CF    UMFCRM(Otl)     VS    UMFCRV<0tl> 

•FEARSCN'     IS    THE    PEARSON    DISTRIBUTION 
BETA    1    =    l.OCt     BETA    2    =    3.80 


NL^eEF    CF    HYPCTHESIS    REJECTIONS     IN    50CO    REPLICATICNS 
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CHI     SQUARED    tpst 
20    PTS  30    PTS 


50    PTS 


CCNTFGL 

10* 
5? 
IS 

371 
206 

27 

431 
237 

48 

464 
227 

45 

FEAFSCN 
VS 

r*  (CD 

10* 
54 
1? 

2231 

1745 

763 

3C39 
2282 
1242 

4476 
3936 
2652 

CCNTFOl 


CCIVTFCL 


1CS 
5? 
IS 


FEARSCN       10? 

h   <c,n       l? 


10? 
5? 
1? 


FE4PSCN       1C3 

VS  53 

N    <  C  ,1 >  I? 


KOLMCG-SMIR    TEST 
20    PTS  3C    °TS 


498 

250 

42 

48  2 

261 

48 

1153 
753 
231 

1481 

1005 

353 

FOUTZ  FN 

TEST 

20  PTS 

30  °TS 

490 

258 

44 

498 

255 

45 

1872 

1220 

469 

2587 

1920 

886 

50  PTS 

517 

261 

54 


2C82 

1447 

621 


50  PTS 

466 

248 

61 


3548 
2894 
1644 


•CCNTFCL'  IS  THE  TEST  CF  UNIFORM<0,1>  VS  UNIFCRMOtl) 

•FEARSCN*  IS  THE  PEAFSON  DISTRIBUTION 
eETA  1  =  2.0C,  BETA  2  =  5.50 


MMEEP  OF  HYPOTHESIS  REJECTIONS  IN  5000  REPLICATIONS 
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